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ABSTRACT 

The  fundamental  concept  of  orthogonality  of  mathematical  objects  occurs  in  a  wide 
variety  of  physical  and  engineering  disciplines.  The  theory  of  orthogonal  functions,  for 
example,  is  central  to  the  development  of  Fourier  series  and  wavelets,  essential  for  signal 
processing.  In  particular,  various  families  of  classical  orthogonal  polynomials  have  tra- 
ditionally been  applied  to  fields  such  as  electrostatics,  numerical  analysis,  and  many 
others. 

This  thesis  develops  the  main  ideas  necessary  for  understanding  the  classical  theory 
of  orthogonal  polynomials.  Special  emphasis  is  given  to  the  Jacobi  polynomials  and  to 
certain  important  subclasses  and  generalizations,  some  recently  discovered.  Using  the 
theory  of  hypergeometric  power  series  and  their  ^—extensions,  various  structural  prop- 
erties and  relations  between  these  classes  arc  systematically  investigated.  Recently,  these 
classes  have  found  significant  applications  in  coding  theory  and  the  study  of  angular 
momentum,  and  hold  much  promise  for  future  applications. 
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I.     INTRODUCTION 

The  abstract  concept  of  orthogonality  of  functions  (or  other  mathematical  objects) 
is  a  generalization  of  the  notion  of  having  two  or  more  vectors  perpendicular  to  one 
another.  This  concept  arises  naturally  in  a  wide  variety  of  physical  and  engineering 
disciplines.  For  example,  the  theory  of  orthogonal  functions  is  central  to  the  develop- 
ment of  Fourier  series  and  wavelets  which  are  essential  to  signal  processing. 

Classical  Fourier  series  (real  form)  depend  on  the  property  that  the  trigonometric 
functions  sine  and  cosine  are  orthogonal  (on  an  appropriate  real  interval)  in  a  formal 
sense  that  will  be  made  precise  later.  As  a  consequence,  a  bounded  periodic  function 
f(x)  of  period  In  which  satisfies  the  Dirichlet  conditions!  may  be  expressed  in  the  form 


f{x)  ~  —  a0  +  y    (  an  cos  nx  +  bn  sin  nx  ) 


/?=i 


where 


1     f* 
an  =  —        f(x)  cos  nx  dx  ,    n  =  0,  1 ,  2, 

bn  =  —        f{x)  sin  nx  dx  ,    n  =  1, 2,  3, . 


are  the  classical  Fourier  coefficients.  These  formulas  can  be  modified  via  a  change  of 
variable  to  accomodate  any  such  function  of  period  2L.  [Ref.  1:  p.  529] 

This  property  can  be  used  to  generate  other  classes  of  orthogonal  functions  - 
polynomials,  for  example  -  that  behave  in  very  structured  and  useful  ways  such  as  in 
generalized  Fourier  series.  In  particular,  specific  families  of  these  "classical  orthogonal 
polynomials"  have  traditionally  been  used  for  solving  problems  arising  in  various  areas 
of  applied  mathematics,  physics,  and  engineering,  among  others. 

This  thesis  develops  the  main  ideas  necessary  for  understanding  the  classical  theory 
of  orthogonal  polynomials.  Special  emphasis  is  given  to  the  Jacobi  polynomials  and  to 
certain  important  subclasses  and  generalizations.     Much  of  the  investigation  will  be 


1  Dirchlet  conditions:   (i)  In  any  period /(.t)  is  continuous,  except  possibly  for  a  finite  number 
of  jump  discontinuities,  (ii)  In  any  period  f{x)  has  only  a  finite  number  of  maxima  and  minim. i. 


made  using  the  theory  of  hypcrgeometric  power  series  and  their  q  —extensions.  The 
classes  discussed  in  this  thesis  are  but  a  small  fraction  of  those  identified  and  studied  in 
the  literature. 

A.     CHEBYSHEV  POLYNOMIALS 

The  Chebyshev  polynomials  of  the  first  kind,  Tn(x),  arise  from  an  elementary  trigono- 
metric consideration.  As  such,  they  satisfy  various  properties  and  identities  which  are 
easily  derived  directly  from  their  definition,  many  of  which  are  observable  from  their 
graphs  (see  below).  This  class  of  polynomials  will  serve  as  the  model  for  some  of  the 
basic  structure  of  more  general  classes. 

For  n  =  0,1,2, ...  ,  define 

Tn(x)  =  cos(a7  arccos  x),    —  1  <,  x  <,  1 

i.e.,  letting  x  =  cos  0  ,  0  <;  6  <  n, 

(1)  Tn(  cos  0)  =  cos  nO. 

Some  immediate  consequences  of  (1)  are  |r„(jc)|<l  for  | jc |  <  1,  with 
T„( cos-—")  =  (-1)*.  ^<k<n;  in  particular  Tn(l)  =  1  and  Tn{ -1)  =  ( -1)"  for  all  n 
which  can  be  seen  graphically  in  Figure  1. 

1.     Three-term  Recurrence  Relation/Differential  Equation 

From  (1),  we  have 

(2)  7o(x)=l  and  T,{x)  =  x, 

and  by  considering  the  identity 

(3)  cos(a  +  b)  +  cos(a  —  b)  =  2  cos  a  cos  b 
with  a  =  nO,  b  =  0,  we  obtain 

(4)  rn+1(*)-2*7;(*)  -  7;_,(.v). 

Equation  (4)  is  known  as  the  three-term  recurrence  relation  for  Tn(x)  which  to- 
gether with  initial  conditions  (2)  imply  that  Tn{x)  is  a  polynomial  of  degree  exactly  n, 
called  the  n"1  Chebyshev  polynomial  of  the  first  kind.  Note  that  the  leading  coefficient 
of  T„(x)  is  2n_l  for  n  >  2.  An  inductive  argument  applied  to  this  recursion  shows  that 


Tn(x) 


-i 


Figure  1.      Chebyshev  Polynomials 

T„(x)  is  an  even  function  if /t  is  even,  and  odd  if  n  is  odd  (see  Figure  1).  The  first  few  are 
listed  below: 

r0(*)-l,    Tl(x)=x,    T2(x)  =  2x2-l,    T3(x)  =  4x2  -  3* 

r4(jc)  =  8;c4  -  8;t2  +  1  ,    T5(x)  =  16jcs  -  20jc3  +  5jc. 

Differentiating  (1)  twice  with  respect  to  &  yields  the  second  order  differential 
equation   for  Tn{x)  : 


(5) 


(l-x2)7V'W  -  xTn'(x)  +  «2r„(jf)-0. 


2.     Orthogonality  of  Chebyshev  Polynomials 

Let 


,n  =       cos 


Imr=  I    cos  mO  cos  nOdB. 


If  m  #  «,  then  using  (3)  with  a  =  m6,  b  =  nd  yields 


I     "  sin(m  +  n)0  sin(w  —  n)6 


m  —  n 


lm/i       2  |         m  +  n         + 
If  m  =  n  i=-  0,  then  by  using  the  identity  cos2a  =    -7  (1  +  cos  2a), 


=  0. 


1  =-r 

"•"    2  L 


d  +   Tr-sin2 

2/7 


If  m  =  n  =  0,  then 


Hence, 


),o  =      <#  =  *■ 

Jo 


(6) 


cos  mO  cos  nd  dd  =  hn    5mn 


where 


-1 


/c  - 


n,      m  =  n  =  0 
tt/2,    m  =  ni^Q 


and 


J0' 


mi=  n 


in  =  // 


is  the  Kronecker  delta  function.  Changing  variables  via  jc  =  cos  0,  we  have 


(7)  f '  Tm(x)Tn(x)      j=L=  dx  =  h~'  5m>n. 

-i  v1_-v 

This  important  property  is  formally  known  as  the  orthogonality  relation  for  the 
Chebyshev  polynomials.  The  reason  for  this  terminology  will  become  clear  in  the  next 
chapter. 

3.  Zeros  of  Chebyshev  Polynomials 

2k  —  1 
Setting       T„(  cos  6)  =  cos  nd  =  0,      we       obtain      d  =  dkn  = — n,      i.e., 

x  =  xk>n  =  cosOk<n,\  <k<n. 

Thus  all  the  zeros  of  Tn(x)  are  real,  distinct,  and  may  be  regarded  as  the 
projections  onto  the  interval  (—1,1)  of  the  equally  distributed  points  Bkn  on  the  unit 
circle,  as  seen  in  Figure  2.  Moreover,  for  1  <  k  <  n,  an  easy  algebraic  check  verifies  that 
0*,  n+i  <  #*,  n  <  0*+i,  «+i  and  therefore  xk<n+x  <  xk<  „  <  xk+Un+l.  Hence,  the  zeros  of  Tn+l(x) 
interlace  with  those  of  Tn(x).  This  interlacing  of  zeros  is  a  striking  feature  of  the  plots 
in  Figure  1. 

The  zeros  of  Chebyshev  polynomials,  and  of  other  orthogonal  polynomials  in 
general,  are  extremely  important  for  applications  to  numerical  analysis,  electrostatics, 
and  many  other  fields. 

4.  Looking  Ahead 

Many  of  the  properties  derived  for  the  Chebyshev  polynomials  Tn(x)  from  their 
trigonometric  definition  (1),  extend  to  more  general  classes  of  orthogonal  polynomials 
via  a  general  theory,  elements  of  which  will  be  developed  in  this  thesis.  Some  of  the 
many  properties  satisfied  by  these  classes  that  we  will  derive  include: 

1.  Orthogonality  with  respect  to  a  weight  function 

2.  Three-term  recurrence  relation 

3.  Second  order  differential  or  difference  equation 

4.  Hypergeometric  series  expression 

5.  Rodrigues'  formula 

6.  Generating  function. 

The  general  approach  we  will  take  is  to  define  these  "classical  orthogonal  polynomials" 
via  terminating  hypergeometric  power  series,  and  from  this  prove  (most  of)  the  other 
properties.  However,  because  of  this  equivalence,  many  authors  choose  to  define  a  given 
class  using  one  of  these  other  characterizing  properties. 


Figure  2.      Zeros  of  Chebyshev  Polynomials 


In  order  to  understand  the  general  theory,  it  is  first  necessary  to  define  the  ab- 
stract concept  of  orthogonality  in  an  appropriately  defined  "space"  of  functions.  We 
turn  our  attention  to  these  fundamental  ideas  in  the  next  chapter. 


II.     BACKGROUND 

A.     ELEMENTARY  LINEAR  ALGEBRA 
1.     Vector  Spaces 

Let  W  denote  the  collection  of  all  vectors  (n-tuplcs),  u  =  (au  a2, ...  ,  a„),  where 
each  ai  e  R,  i  =  1,2, ...  ,  n.  The  standard  Euclidean  inner  product  (also  referred  to  as  the 
dot  product  )  of  two  such  vectors  u  =  (au  a2, ...  ,  a„)  and  v  =  (bu  b2, ...  ,  bn)  is  given  by 


n 

<u,  v>  =2^ci<b 


i- 


The  length,  or  norm,   of  a  vector  u  e  R"  is  given  by 

||u||  =  7<u,u>    =ma,l 

Two  vectors  u,  v  e  R"  are  perpendicular,  or  orthogonal,   if  and  only  if  <  u,  v>  =  0. 

The  objective  of  this  chapter  is  to  extend  these  familiar  notions  to  objects  other 
than  classical  Euclidean  vectors,  in  particular,  the  "vector  space"  of  polynomials  defined 
on  a  real  interval  [a,  b~\. 

A  vector  space  V  over  a  scalar  field  F  (usually  R  or  C  )  is  a  nonempty  set  of 
objects  called  vectors,  for  which  the  operations  of  addition  and  scalar  multiplication  are 
defined.  Addition  is  a  rule  for  associating  with  each  pair  of  vectors  u  and  v  in  V  an 
element  u  +  v,  called  the  sum  of  u  and  v  .  Scalar  multiplication  is  a  rule  for  associating 
with  each  scalar  c  in  F  and  each  vector  u  in  V  an  element  cu,  called  the  scalar  multiple 
ofu   bye.   [Ref.  2:  p.  150] 

For  all  u,  v,  vv  e  V  and  c,  d  e  F,  a  vector  space  V  must  satisfy: 

1.  Additive  closure.   u,\eV=>u  +  \eV 

2.  Commutativity.   u  +  v  =  v  +  u 

3.  Associativity,   u  +  (v  +  w)  =  (u  +  v)  +  w 

4.  Additive  identity.   There  exists  a  zero  vector,  0  e  V,  such  that  0  +  u  =  u  4-  0  =  u. 

5.  Additive   inverse.      For   each   u  e  V,   there   exists   a   vector  —  u  e  V,   such   that 
u  +  (  -  u)  =  (  -  n)  +  u  =  0. 

6.  Multiplicative  closure,    u  e  V  and  c  e  I;  =>  cu  e  V 


7.  Distributivity.   c(u  +  v)  =  cu  +  cv 

8.  Distributivity.   (c  +  d)u  =  cu  +  du 

9.  Multiplicative  associativity.   c(</u)  =  (cd)u 

10.  Multiplicative  identity.   There  exists  a  scalar  1  e  F  such  that  lu  =  u. 

Example  1:    W  (the  model) 

Example  2:    PN\_a,b~\  =  {polynomials  of  degree  <  N  on  the  interval  [a,/;]} 

Example  3:    P[a,b~\  =  {polynomials  on  the  interval  [a,/;]} 

Example  4:    C\_a,f\  =  {continuous  functions  on  the  interval  [a,/']} 

Note  that  PN\_a,b~\  c  P[a,b"\  <=  C[«,&].  These  vector  spaces  are  sometimes  re- 
ferred to  as  function  spaces.  The  interval  \_a,b~\  may  be  finite  or  infinite  (i.  e., 
\_a,  oo) ,  (— oo,  6]  ,  or  (— oo,  oo)  )  for  our  purposes. 

A  subset  U  of  a  vector  space  V  is  said  to  be  a  vector  subspace  of  V  if  it  is  a 
vector  space  in  its  own  right. 

Example  5:  PA[a,b~]  is  a  subspace  of  P[_a,b~],  which  in  turn  is  a  subspace  of 
Cla,bl 

Given  a  set  of  vectors  {v,,  v2, ...  ,  v„}  in  a  vector  space  V,  and  scalars 
c,,  c2, ...  ,  cn,  the  vector  c1v1  +  c2\2  +  •••  +  cn\n  e  V  is  said  to  be  a  linear  combination  of 
{v,,  v2f ...  ,  vj.  The  set  of  vectors  {vM  v2,  ...  ,  v„}  is  said  to  be  linearly  dependent  if  there 
exist  scalars  cu  c2, ...  ,  c„,  not  all  equal  to  zero,  such  that  the  linear  combination 
civi  +  C2V2  +  '"  +  cnsn  —  0-  ( Equivalcntly ,  at  least  one  of  the  vectors  v,  can  be  expressed 
as  a  linear  combination  of  the  others.)  Otherwise,  the  set  {v„  v2,  ...  ,  vj  is  linearly  inde- 
pendent. An  infinite  set  5  =  {v,,  v2, ...  ,  v„, ...  }  is  defined  to  be  linearly  independent  if  every 
finite  subset  of  S  is  linearly  independent;  otherwise  S  is  linearly  dependent  [Ref.  3:  p.  8]. 
The  vectors  {v,,  v2, ...  ,  v„}  are  said  to  span  V  if  every  vector  \  e  V  can  be  represented  as 
a  linear  combination  of  {v„  v2, ...  ,  v„}.  In  this  case,  we  write  V '=  span{\x,  v2, ...  ,  v„}.  The 
vectors  {v,,  v2, ...  ,  v„}  form  a  basis  for  V  if  they  are  linearly  independent  and  span  V.  The 
dimension   of  V  is  the  number  of  elements  in  any  basis. 

Example  6:  The  set  {e,,  e2,  ...  ,  e„}  is  the  standard  basis  for  IR",  where 
e,  =  (0,  0, ...  ,0,  1,0, ...  ,  0)  i.e.,  the  vector  with  a  one  in  the  i'h  position  and  zeros  else- 
where, /=  1,  2,  ...  ,  n. 

Example  7:  The  set  {1,  x,  x2,  ...  ,  .rv}  is  the  standard  basis  for  PJatb'].  (Linear 
independence  is  ensured  by  the  r:undamental  Theorem  of  Algebra.)  The  dimension  of 
PN[a,b]  is  therefore  N+l. 


Example  8:    The  set  {1,  x,  x2, ...  ,  x", ...  }  is  the  standard  basis  for  P[a,6],  and 
hence  P\_a,b~\  is  an  infinite-dimensional  vector  space. 
2.     Inner  Product  Spaces 

An  inner  product  on  a  real  vector  space  V  is  a  mapping 

<  ,   >:Vx  V^U 

such  that  for  all  u,  v,  w  e  V  and  a,  /?  e  U,  the  following  properties  hold: 

1.  Positive  definiteness:    <  u,  u>  >  0,  and  <  u,  u>  =  0  if  and  only  if  u  =  0 

2.  Symmetry:    <  u,  v  >  =  <  v,  u  > 

3.  Bilinearity:    <au  +  /?v,  \v>  =  a  <  u,  \v>  +  /?<v,  w> 

A  vector  space  with  an  inner  product  is  known  as  an  inner  product  space. 

Example  9:    V  —  U"  ;  let  constant  "weights"  w,  >  0  be  given,  /=  1,  2, ...  ,  n. 
For  u  =  (als  a2, ...  ,  an)  and  v  =  (blt  b2, ...  ,  b„),    u,  v  e  V, 

n 

<  u  ,  v  >  =  y  at  b{  wt 

If  wt  =  1  for  /=  1,  2, ...  ,n,  then  this  reduces  to  the  standard  Euclidean  inner  product, 
or  dot  product.   Otherwise,  this  is  referred  to  as  a  weighted  inner  product. 

The  next  two  examples  are  commonly  applied  inner  products  on  function  space, 
and  are  analogues  of  the  previous  example.     We  assume  a  given  weight  function 
w(x)  >  0  in  (a, b),  integrable  in  the  first  case  (e.g.,  continuous  for  [a,b~\  a  finite  interval). 

Example  10:    V=  PNla,b~\,  P[a,b~],  C{a,b] 

<f,g>  =  \   f(x)g(x)w(x)dx 

Example  11:    V  =  PN{.a,b1 

N 

</,g>  ^Yj/wsixx*) 

(Positive  definiteness  is  ensured  by  the  Fundamental  Theorem  of  Algebra.) 


The  norm  induced  by  the  inner  product  is  given  by  ||u||  =  J  <  u,  u>  .2 
Example  12:    For  the  inner  products  of  Examples  10  and  11  therefore 

11/11  =  (  j\f(x)]2w(x)dx  V2      and 


V  ^ 

1/2 


11/11=       ^t/tof"**) 


\  x=0  / 

respectively.   These  are  sometimes  referred  to  as  "L2-norms." 

Two  vectors  u,  v  e  V  are  said  to  be  orthogonal,  denoted  ulv,  if  and  only  if 
<  u,  v>  =0.  The  vectors  u  and  v  are  said  to  be  orthonormal  if  ulv  and 
IMI  =  IMI  =  1-  Note  that  the  orthogonality  of  vectors  in  a  space  is  determined  by  the  in- 
ner product  being  used. 

The  two  examples  which  follow  refer  back  to  Chapter  I,  Section  A.l. 

Example  13:  Formula  (6)  shows  that  the  functions  {  1,  cos  x,  cos  2x, ...  }  are 
orthogonal  on  [0,  n]  with  respect  to  the  uniform  weight  function  w{x)  =  1.  A  similar 
computation  shows  that  the  same  property  holds  on  [— n,  n~]  with  respect  to  the  weight 
function  w(x)  =  —  ,  i.e., 

<f,g>  —jrj*f(x)g(x)dx. 

One  advantage  of  preferring  this  inner  product  over  the  standard  one  lies  in  the  com- 
putation of  norms.  Using  ||/||  =  J  <  f,  f  >  ,  we  have  ||  1  ||  =  2  and  ||  cos  nx  ||  =  1  if 
n  >  1.  Hence  the  functions  {  1/2,  cosjc,  cos  2x, ...  }  are  orthonormal  on  [—re,  n]  with 
respect  to  the  inner  product  above.  Similar  statements  hold  for  the  integral  of  a  product 
of  two  sine  functions  on  [—re,  it],  as  well  as  for  the  product  of  a  sine  and  a  cosine. 

Example  14:  By  (7)  the  Chebyshev  polynomials  {T„(x)}  form  an  orthogonal 
class  with  respect  to  the  inner  product  of  Example  10  above  on  [  — 1,1]  with  the  weight 
function  w(x)  =  (1  —  x2)'1'2. 


2  Recall  that  <  u,  u  >  >  0.   We  remark  that  in  the  same  way  we  defined  inner  product  earlier, 
it  is  possible  to  define  a  general  norm  on  a  vector  space  which  is  not  induced  by  an  inner  product. 
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We  remark  here  that  Examples  10  and  11  can  be  unified  into  a  single  inner 
product  on  a  "polynomial"  space  V  via 


/,*>  =  \f(x)g(x)da(x), 


< 

where  da.{x)  is  a  positive  Lebesgue-Stieltjes  measure  on  a  measurable  set  E  possessing 
finite  moments,  i.e.,  x"  da(x)  integrable,  n  =  0,  1,2,....  In  Example  10,  £=  [_a,b]  a  R  and 
doc(x)  =  w{x)dx\  the  resulting  expression  is  known  as  a  continuous  inner  product,  while 
in  Example  11  the  set  E  consists  of  a  finite  number  of  points  {0,  1, ...  ,  iV}  elR,  and  the 
associated  measure  gives  rise  to  a  discrete   inner  product. 

B.     FOURIER  SERIES 

Let  v  e  V,  and  U  be  an  ^-dimensional  subspace  of  Shaving  some  orthonormal  basis 
{  Ui, ...  ,  u„ }.  (Any  basis  can  be  orthonormalized  via  the  Gram-Schmidt  process  -  see  next 
section.)  The  vector  v  can  be  resolved  into  a  sum  of  two  vector  components: 

(1)  v  =  (v  —  w)  -f  w 

where  wet/  and  (v  —  >v)  1  U.  (See  Figure  3.)  The  vector  w  is  referred  to  as  the 
orthogonal  projection  of  v  onto  U.  Since  the  vector  (v  —  w)  is  orthogonal  to  every  vector 
in  U  by  construction,  it  follows  that  for  eachy=  1,  2, ...  ,  n,  <  v— w  ,  u,  >  =  0,  or 


(2)  <  v  ,  u;  >  =  <  w ,  u; 


> 


Moreover,  since  it  lies  in  U,    vector  w  can  be  expressed  as  some  linear  combination  of 
{  u„  ... ,  u, }  : 


n 


Take  the  inner  product  of  both  sides  with  u,  for  each  j  =  1,  2, ... ,  n.  From  the  assumption 
that  <  u, ,  u,  >  =  0  unless  i=j,  we  have  the  property  that 

<   W  ,    Uj  >   =Cj   <   Uj,    Uj>  . 

Thus, 

(3)  cj  =  <  v  ,  uj  > 
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via  (2)  and  the  assumption  that  <  u, ,  u,  >  =  ||  u,  ||2  =  1. 
Thus, 

n 

(4)  w  =  2_j  K  v  '  ui  >  u< 

and  this  vector  represents  the  "best  approximation"  in  U  to  v  e  V  in  the  sense  that  of 
all  vectors  z  e  U,  it  is  the  projection  vector  >v  e  U  which  uniquely  minimizes  the  distance 
l|v-z||. 

Suppose  now  that   U     is  an  infinite-dimensional  subspace  of  V     (also  infinite- 
dimensional),  having  orthonormal  basis  {  u,, ...  ,  u„, ...  }.  Then  from  (1)  and  (4) 


=  (v-vv)  +  2^  <  v  ,  u,  >  u. 


/=1 


we  may  write 


oo 


(5)  v  =  2,  <  v  ,  u,  >  U/ 

(=1 

in  the  sense  that 

(6)  lim  ||  v-w  ||  =  0, 

i.e.,  the  norm  of  the  "residual  vector"  (and  hence  the  vector  itself)  v— w  -*  0  as 
n  -*  oo.  Formula  (5)  is  known  as  the  generalized  Fourier  series  for  v  e  V  with  respect 
to  the  orthonormal  basis  {  u, }°!,.  The  coefficients  given  in  (3)  are  called  the  generalized 
Fourier  coefficients  of  ve  V.  Statement  (6)  is  known  as  the  norm  convergence  property 
of  Fourier  series,  and  the  "minimization  property"  mentioned  above  extends  to  this 
infinite-dimensional  case. 

Example  15:  Let  V=C[a,b],  and  { <t>,(x)  }°T„  be  an  orthonormal  basis  of 
eigenfunctions  (sometimes  referred  to  simply  as  an  eigenbasis  )  of  V.  Then/  e  V  has  a 
Fourier  scries  representation 


i=0 
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Figure  3.      Best  Approximation 


with  Fourier  coefficients 


(7) 


:,  -  <  /.  0i  >  -  [  /(*)  U*)  WW  dx. 

•J  SI 


In  this  function  space  context,  norm  convergence 


lim 

n— »oo 


i=0 


2  \l/2 

dx  )      =0 
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is  referred  to  as  mean  square  convergence,   and  is  the  least  squares  principle  in  regression 
analysis. 

In   particular,   if   V     is   equipped  with   the   inner   product   of  Example    13   and 
orthonormal  basis 

{  0/jc)  }  =  {  1/2,    cos  x,    cos  2x, ...  ,    sin  x,    sin  2x, ... } 

on  [— n,  7r]  (see  Example  13,  Section  A. 2),  then  a  suitable  function  fe  V  (and  its 
2^-periodic  extension  on  IR  )  has  a  classical  Fourier  series 


OO 

f(x)  ~  —  a0  +  /.(an   cos  nx  +  °n   sm  nx  ) 

n=0 


where 


1     f* 

an  =  <  f,    cos  nx  >  =  —        f(x)  cos  /7jc  <^c 

1     f* 
bn—  <  f,    s'mnx  >  =—        f(x)  sin  /7„v  flte, 

as  indicated  in  the  Introduction. 

C.     GRAM-SCHMIDT  ORTHONORMALIZATION 

The  Gram-Schmidt  process  orthonormalizes  any  set  of  linearly  independent  vectors 
in  an  inner  product  space.  This  method  will  be  used  in  later  sections  for  different  inner 
products  on  the  vector  space  P\_a,b~]. 

Begin  with  an  inner  product  space  V  and  any  set  of  vectors  {v,,  v2, ...  ,  v„, ... },  finite 
or  infinite,  such  that  any  finite  number  of  elements  of  this  set  are  linearly  independent. 
Recursively  define  a  new  set  of  vectors  {u,,  u2, ...  ,  u„, ...  } 


Uk  =  tj'    k=^2—^ 


where  yk  =  yk-  \\k,  vvith3 


3  By  convention,  Va,  =  0,  giving  w,  =  0. 
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k-\ 
>V^=2j<Vfc,Uf>U/. 


l=\ 


These  new  vectors  {u,,  ti2, ...  ,  un, ...  }  are  orthonormal  by  construction  and  span  the 
same  space  as  the  original  vectors.  Note  that  this  process  occurs  in  two  stages: 
orthogonalization  and  normalization.  The  orthogonalization  is  accomplished  by  sub- 
tracting y/k,  the  orthogonal  projection  of  \k  onto  the  subspace  spanned  by 
{u„  u2, ... ,  u*_j}.  The  component  of  \k  which  remains,  denoted  above  as  yk,  is  then 
orthogonal  to  the  vectors  {u,,  u2, ...  ,  u^_,}  as  shown  in  Figure  4.  The  normalization  is 
then  achieved  by  dividing  yk  by  its  norm,  thus  giving  it  unit  "length". 
1.     Legendre  Polynomials 

Example  16:  Let  V=P\_— 1,1]  with  basis  {l,.v, x2, ... ,  x", ...  }  and  uniform 
weight  function  w(x)  =  1.  The  inner  product  is  then  given  by  </, g  >  =  /'  f(x)g(x)dx. 
The  Gram-Schmidt  process  yields  the  set 


{4=i  =  {^wu=  -f-vy^  ~~ r~"vv  _TJ'  "^r~Vv  ~y-V"-( 

as  an  orthonormal  basis  for  P[— 1,1].  Since  this  set  is  linearly  independent,  we  can 
standardize  the  set  by  taking  scalar  multiples  of  these  polynomials  so  that  Pn{\)=  1. 
Members  of  the  resulting  orthogonal  set 

{Pn(x))Zo  =  {l^-y  0*2  -  1), y  (5-v3  -  3*),  ...  } 

are  known  as  the  Legendre  polynomials  on  [—1,1].  If  the  normalized  Legendre 
polynomials  { pn{x)  }^10  are  used  as  the  orthogonal  eigenbasis  for  a  Fourier  series,  the  re- 
sulting expansion  is  often  referred  to  as  a  Legendre  series  representation;  when 
Chebyshev  polynomials  are  used,  we  obtain  a  Chebyshev  series  representation,  etc. 

The  Gram-Schmidt  process  can  always  be  used  in  this  way  to  generate  a  class 
of  orthogonal  polynomials  with  respect  to  a  given  inner  product  (i.e.,  weight  function) 
on  a  real  interval.  When  using  the  Gram-Schmidt  process  from  the  basis 
{l,x,  x2, ...  ,xn, ...  },  the  orthogonalization  stage  producing  y*  results  in  a  set  of  monic 
polynomials,  i.e.,  the  leading  coefficient  of  each  polynomial  is  one.  In  the  normalization 
stage,  we  are  dividing  by  the  norm  ||  yk  ||  >  0.  Thus  the  leading  coefficient  of  polynomials 
in  an  orthogonal  class  is  strictly 'positive.     In  the  next  chapter,  we  will  examine  other 
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Figure  4. 


Orthonormalization 


ways  to  define  these  classes.     It  is  the  structure  and  applications  of  certain  of  these 
classes  with  which  we  will  primarily  be  concerned. 

D.    THE  GAMMA  FUNCTION 

The  gamma  function  V{x)  is  a  fundamental  mathematical  object  that  appears  fre- 
quently in  the  representations  of  orthogonal  polynomials  as  well  as  in  many  other  ap- 
plications. This  "special  function"  was  developed  as  a  generalization  of  the  factorial 
function  of  the  natural  numbers.  As  we  will  see,  the  gamma  function  has  the  value 
(n  —  1)!  for  the  positive  integers  n  but  it  is  defined  for  noninteger  values  as  well. 

A  conventional  definition  for  the  gamma  function  is 
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(8) 


r(x)  = 


e~'tx~xdt,  x>0. 


The  positivity  of  x  ensures  that  this  improper  integral  converges.  We  now  develop  some 
fundamental  properties  of  the  gamma  function.    Integration  by  parts  in  (8)  yields 


(9) 


r^+i)  =  ^r(.v). 


We  now  introduce  the  Pochhammer  symbol  or  shifted  factorial,  {a)n ,  to  simplify  our 
notation.    For  n  >  0,  define 

(a)n  =  a{a  +  l)(a  +  2) ...  (a  +  n  -  1)  ,    if  n  >  1 

and  {a)0  =  1.  Letting  a  =  1  gives  (1)„  =  (1)(2)(3) ...  (n)  =  n\.  Note  that  for  a  negative  inte- 
ger, (— m)n  =  0  if  n  >  m  >  0.  The  shifted  factorial  can  be  defined  for  negative  subscripts 
but  we  will  not  need  this  in  our  work  with  polynomials. 
Iteration  of  (9)  n   times  yields 


(10) 


T{x  +  n)  =  {x)nr{x) 


for  every  positive  integer  n.    Using  this  property,  the  gamma  function  can  be  extended 
to  include  negative  real  numbers  by  defining 


(11) 


r(.v)  =  — —  r{x  +  n)    for    —  n  <  x  <  —  n  +  1 . 


Since  this  expression  is  undefined  when  x  is  zero  or  a  negative  integer,  the  gamma 
function  is  not  defined  for  those  values. 

Letting  x  =  1  in  (8)  and  computing  directly,  we  have  T(l)  =  1.  It  then  follows  that 
r(n  +!)  =  «!  by  letting  x  =  1  in  ( 10).    Furthermore, 


f*oo 


-'r{]l2)dt  =  2 


2  I 

e  "  du  =  J  7i 


where  the  second  integral  can  be  evaluated  by  standard  methods  involving  multiple  in- 
tegrals. 
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Finally,  we  define  a  generalized  binomial  coefficient  as  follows.    For  .v    and  a  non- 
negative  integers,  define 


x+a\       (x+a)\ 
x 


For  nonintegral  a,  define 


jc+a 
x 


(a+\)x        ru+a+i) 


(i),       ru+i)r(a+i) 


1.     The  Beta  Function 

An  integral  related  to  the  gamma  function  defines  another  useful  function  called 
the  beta  function  which  is  given  by 


P\ 


(12) 


B(x,y)  = 


tx-\\-t)y-xdt 


Jo 


for  x  ,y>Q.  We  now  establish  an  important  connection  between  the  beta  and  gamma 
functions.   We  start  with  an  identity  easily  verified  from  (12): 


(13) 


Also  from  (12) 


B(x,  y+l)  =  B(x,  y)-B(x+i,  y). 


B(x+l,y)  = 


t\\-t)y-xdt 


which  when  integrated  by  parts  gives 


x 


B(x+l,y)  =  ^B(x,y+l) 


y 


Substituting  into  (13),  we  obtain 


B(x,y)  =  ^rLB(x,y  +  1) 


which  when  iterated  vields 
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#(*,  >')  =  —r~\ —  #(-v*  ^  +  n)  =  — r~; —     r     ( J  _  0       dt- 


[y)n 


(y)n 


Changing  variables  from  t  to  —  , 


(14) 


Taking  the  limit  as  n  -*  oo,  and  using  the  fact  that  lim(  1 — —  )"  =  e~',  we  have 

(x+y)n 


(15) 


B(x,  y)  =  r(x)  lim 


^°°  (y)«« 


(The  fact  that  we  can  pass  the  limit  through  the  integral  on  the  right  can  be  mathemat- 
ically justified.)   Ify  =  1,  then  (15)  gives 

B{x,l)  =  f(jc) Jim  J"  . 


By  direct  evaluation  using  (12), 


Hence 


which  can  be  written  as 


n 


D(x,l)  = 


r-'A-i 


•'O 


1       mr       (x+1)* 
—  =  r(.t)Jirn        ,    x     , 


n-*00  yi\ 


n:  n 


T{x)  =  lim 


n\  n 


n~*°°  x(x+  1)„  ' 


Noting  that  x(x  +  l)„  =  (x)n^  =  {x)n(x  +  n),  we  have 
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x-\ 


n\  n      n 


Y{x)  =  lim 

«-°o  (x)n(x  +  n) 


which  gives  the  form 

n\  n 


(16)  r(.v)  =  lim 

«-oo        (*)„ 

(Equation  (16)  was  Euler's  original  definition  of  the  gamma  function.  A  separate  "esti- 
mation" argument  may  be  used  to  show  that  this  limit  mathematically  exists.)  Thus  by 
(15), 

!     x+y-\ 

D(x,y)  =  r(x)\[m-^- . 

I     7-1 

Then  by  (16),  we  have  the  useful  identity 

(17)  IK*,  j-)-  r{x+y)  ■ 

We  will  find  this  identity  useful  in  understanding  the  Jacobi  polynomials  in 
Chapter  IV,  where  it  becomes  necessary  to  evaluate  a  related  integral: 


I     [1-xf  {1+xf  dx. 


We  remark  here  for  future  reference  that  the  formal  change  of  variable  x  =  1—  It  can  be 
used  to  transform  this  integral  into 

Jo 

=  a+/3+[  n«+i)r(/?+i) 

r(a+/J+2)       ' 
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III.     GENERAL  THEORY  OF  CLASSICAL  ORTHOGOiNAL 

POLYNOMIALS 

In  this  chapter,  we  examine  some  of  the  characteristic  properties  associated  with 
classes  of  orthogonal  polynomials.  Some  of  these  properties  provide  alternate  means 
of  defining  a  class.  These  alternate  definitions  often  provide  a  straightforward  way  of 
producing  a  specific  result  that  may  be  very  difficult  to  derive  otherwise. 

Throughout  this  chapter,  we  let  {p„(x)}™0  denote  a  set  of  real  polynomials  with 
pn(x)  of  degree  n,    i.e., 

(1)  pn{x)  =  knxn  +  snxn~]  +■■■     ,  kn>Q. 

Recall  that  these  polynomials  are  said  to  be  orthogonal  on  an  interval  [a, b~]  with  respect 
to  a  continuous  weight  function  w(x)  >  0  on  (a,b)  if 

(2)  <  pm  ,pn  >  =     pm(x)pn(x)w{x)dx  =  h~]  Smn 

where  the  normalization  /?"'  ^  0  is  chosen  to  simplify  the  expression  of  certain  formulas. 
Note  that  since 


rb 
Pn  II2  =  J     [  Pni.XJ\w{x)dx  =  h~\ 


it  follows  that  A„  >  0. 


A.     POLYNOMIAL  EXPANSIONS 

We  begin  by  showing  that  any  real  polynomial  qm{x)  of  degree  m    on  [a,//j  can  be 
written  as  a  linear  combination  of  orthogonal  polynomials  {pn(x)}^0  : 

m 

(3)  q„{x)  =  2^  Hm  AW 

t-Q 

for  constants  a.,m  ,  /'  =  0,  1,  2, ...  ,  m. 

The  proof  is  by  induction  on  the  degree  in.      Since  qm{x)  's  a  rcal  polynomial,  we 

write 
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<lm\x)  =  amX      +    bmx  +     - 

where  am  ^  0.  For  m  =  0,  (3)  reduces  to  a0  =  a00  k„  using  the  form  for  pn(x)  given  in  (1) 
and  aOi0  is  uniquely  determined. 

For  the  induction  hypothesis,  suppose  that  for  in  >  1,  we  can  write  any  polynomial 
of  degree  m  —  1  as  a  linear  combination  of  { p„(x)  }™~q  : 

m-\ 
1=0 

Since  qm{x)  —  (ajkm)pm(x)  =  q„.x{x)  is  a  polynomial  of  degree  m  —  1  the  induction  hy- 
pothesis implies  there  is  a  representation 

m-\ 


qm{x)  ~{~j~)  p*W  =  Aj  °Wi  ^W 


Now  set  am>m  =  (ajkm)  and  the  result  (3)  follows.   [Ref.  4:  p.  33] 

Using  the  theory  of  Fourier  series  developed  earlier,  we  next  determine  the  coeffi- 
cients a,m  explicitly.    For  /  =  0,  1, ...  ,  m,  let  c,m  =  a^jjl^  and  let 

Hx)=i^m=y/  iPi{x)- 

Then  by  construction,  { <t>,(x)  }~0  is  an  orthonormal  set  of  polynomials.   Writing  (3)  as 

m 

i=0 

we  see  that  the  right-hand  side  can  be  interpreted  as  a  (terminating)  Fourier  expansion 
of  qjx).  Hence  the  results  of  Example  15  in  Chapter  II,  Section  B  may  be  applied.  In 
particular,  by  (7)  in  that  section,  the  coefficients  c,m  are  given  by 

ci,m  -  <  Urn  >   <f>l  >  =        ?mW  4>fa)  w(x)  dx- 
Changing  back  to  the  old  variables, 
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ai,m  =  hi  <  1m  >  Pi>  =  hi       <7™(-Y)  Pi(x)  w(x)  dx- 
We  are  now  in  a  position  to  show  that  the  orthogonality  property 

(4)  pn(x)pm(x)w(x)dx  =  0    ,m^n 

can  be  expressed  equivalently  as, 


(5) 


b 

pn{x)  xmw(x)dx  =  0  ,  m  <  n. 


To  see  (5),  substitute  the  form  (1)  for  pm{x)  into  (4)  where  m<n.  The  linearity  of  the 
integral  gives  (5).  On  the  other  hand,  since  xm  is  a  polynomial,  we  can  write  xm  as  a 
linear  combination  of  the  orthogonal  polynomials,  so  (5)  gives  (4).  Note  that  (5)  implies 
each  pn(x)  is  orthogonal  to  every  polynomial  of  lower  degree.   [Ref.  4:  pp.  33-34] 

B.     THREE-TERM  RECURRENCE  RELATION 

The  three-term  recurrence  relation   is  a  useful  result  which  holds  for  any  three  con- 
secutive orthogonal  polynomials: 

(6)  pn{x)  =  {Anx+ Bn)pn_x{x)-  Cnpn_2{x)    ,n  =  2,3,4,... 

where  An ,  Bn ,  and  C„  are  constants  given  by 

"-n  ^        t.  ,   (    Sn  sn-\    \        ^  ^n  ^n-2     "n-2 


Kn-\  \   Kn  Kn-\    /  i^n-V         n~l 

The  recurrence  relation  is  valid  for  n  =  1  if  /?_,  =  0  with  C,  arbitrary.  In  this  case,  the 
formula  for  A„  also  holds  for  n=\.  (In  the  contrapositive  form,  this  statement  is  a 
powerful  tool  for  showing  that  a  polynomial  set  is  not  orthogonal.)   [Rcf.  5:  p.  234] 

To  prove  this,  we  begin  by  considering  pn{x)  —  {kjkn_x)  jc/v,(.y),  a  polynomial  of  de- 
gree no  greater  than  (n  —  1).  We  expand  it  in  terms  of  the  orthogonal  polynomials 
{Pi(x)}"~Q  via  the  technique  of  the  previous  section  to  obtain 
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(7) 


Pn(x)  ~  -]f—  xPn-\(x)  =  }^  ai,n-\  Pi(x)' 
"~]  (=0 


The  coefficients  a,„_,  are  determined  by 


(8) 


a,„_i  =  -ht <  xpn_x(x) ,  p^x)  >    ,  i  <  n  -  1 


Vi-1 


Because 

<  x  pn_x (x)  ,  pfac)  >  =  \  x  pn_x{x)  pt(x)  w(x)  dx  =  <  pn_x (x)  ,  x pfa)  >  , 

it  follows  that  for  i  <  n  —  3,  x  p,(x)  is  of  degree  no  greater  than  {n  —  2).  Hence 

< pn-\{x)  >  x Pi(x) >  =  <x pn^(x)  ,  pt{x) >  =  0    ,  /  <  n  -  3, 

since  pn^(x)  is  orthogonal  to  every  polynomial  of  lesser  degree.  Thus  the  constants 
«ia-i»  *a*-n  •••  '  a-3."-i  are  all  zero,  leaving  an_2tn_,  and  a„_lin_,.  With  this  knowledge,  (7) 
becomes 


Setting  An  =  kjkn_x  ,  Bn  =  a„_Iin_, ,  and  Cn  =  -  an_^_,  then  rearranging  terms  gives  (6). 
To  determine  C„  explicitly,  we  write 

C„  =  -  a„_2^_,  =  Anhn_2< pn_x(x) , pn_2{x)> 
from  (8).    Since 


n-\ 


xPn-l(-x)  =  kn-2x  "    +  "' 


Kn-l 


Vi-1 


ln-\ 


[W-'  +  -  ] 


n-2 


Pn-](X)+    >    Pj,n-2Pj(x) 


y=0 


we  can  write 
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A, 


n-2 


Cn  =  hn_2  — </>„_,(*) ,  />„_i(-v)  +  2_,  Pjjn-2  P£x) >  • 

Using   the   properties   of  the   inner  product  and  the   orthogonality  property  of  the 
polynomials,  we  obtain 

*n  ,    v  /    s  _  "n-2       An  ><-nkn_2     nn_2 


Cn  =  hn_2- <Pn-i(x),Pn-l(x)>  = 


ln-\  nn-\      An-\  {k„_x)       "n-\ 

as  given  in  (6).   [Ref.  6:  p.  8] 

Into  equation  (6)  we  substitute  the  expanded  forms  of  the  polynomials 
pn(x) ,  p„-i{x) ,  p„_2(x)  from  (1)  and  the  constants  A„ ,  Cn  from  above.  Equating  coefil- 
cients  of  x"~l  gives  Bn.  Since  kn  >  0  and  hn  >  0,  it  follows  that  A„  >  0  and  Cn  >  0  and  the 
proof  is  complete. 

The  nonnegativity  of  the  constants  A„  and  C„  is  important  for  the  converse  of  the 
result  in  (6).  Favard  showed  that  the  existence  of  a  three-term  recurrence  relation  in  the 
form  of  (6)  implies  that  the  polynomials  of  the  set  are  orthogonal  with  respect  to  some 
weight  function  over  some  interval  using  Stieltjes  integration  [Ref.  7]. 

We  observe  that  in  order  to  generate  the  polynomials  of  an  orthogonal  class  with 
the  recurrence  relation  (6),  we  need  the  sequences  of  constants  A„ ,  Bn ,  and  Cn  together 
with  two  of  three  consecutive  polynomials  in  the  class.  Other  techniques  provide  what 
some  authors  call  a  pure  recurrence  relation  requiring  only  two  of  three  consecutive 
polynomials  to  define  the  class,  because  the  constants  as  functions  of  n  are  contained 
explicitly  in  the  recurrence  relation.  The  recurrence  relation  derived  in  Chapter  I,  Sec- 
tion A.l  for  the  Chebyshev  polynomials  is  an  example. 

C.     CHRISTOFFEL-DARBOUX  FORMULA 

The  ChristolTcl-Darboux  formula  is  an  important  identity  which  can  be  derived  from 
the  three-term  recurrence  relation.   The  identity  is 

v^  hnkn   Pn+1(x)pn(y)  -  pn(x)pn+1(y) 


(9)  LhjPMHy)  =  ~k 


x  —  y 

7=0  "+1  ^ 


To  prove  (9),  note  that  from  (6)  we  have 

pj+i(x)  =  (AJ+xx  +  BJ+l)pj(x)  -  CJ+lpj_{(x) 
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which  we  rearrange  to  give 


&* 


( 1 0)  xpj(x)  =  — pJ+ ,  (x)  -  -j —  pj{x)  +  -j —  pj_ ,  (x) 

Aj+\  nj+\  /lj+\ 


and  similarly 


Bj+\      .   .       Cj+i 


(ii)  ypj(y)  =  —, — Pj+i(y)  -  —A — Pj(y)  +  —A — Pj-i(y)- 

Aj+\  Aj+\  /lj+\ 

These  recurrence  relations  are  valid  fory  =  0  if  we  set  C,  =  0  and  Jp_,(.r)  =  p^{y)  =  0. 

Multiply  (10)  through  by  pj(y)  and  (11)  by  Pj(x),  subtract  the  results,  and  then  mul- 
tiply through  by  h;  to  obtain 

hikj  r  i 

hj  (x  -  y)  pj(x)  pj(  y)  =  - —  [  pj+ ,  (x)pj(  y)  -  pj+ ,  ( y)ppc)\ 

Kj+\ 

+  — i —  [  Pi- 1  (x)pj(  y)  -  Pi- 1  ( y)Pj(x)]  ■ 
Kj 

Summing  overj   from  0  to  n  yields  a  telescoping  series 


v-i                          hn  kn 
(x  -y)}hjPj{x)pj{y)  =  -t [pn+i(x)pn(y)  -  pn(x)pn+](y)\ 

from  which  the  identity  (9)  follows.   [Ref.  4:  p.  39] 

Now  subtract  and  add  the  quantity  p„{x)  p„^(x)  to  the  numerator  of  the  right-hand 
side  of  (9)  and  let  j/  tend  to  x  to  obtain  a  limiting  case  of  the  Christoffel-Darboux  for- 
mula: 

Eh    k. 
hj  [ Pj{x)  f  =  -^-JL  [ pn+]  (x)  pn(x)  -  P'n(x)  pn+l(x)  ]. 

We  will  use  (12)  in  the  next  section. 

D.     ZEROS  OF  ORTHOGONAL  POLYNOMIALS 

In  Chapter  I,  we  observed  that  the  zeros  of  the  Chebyshcv  polynomials  [Tn(x)}  are 
real,  distinct,  and  lie  in  the  interval  (—1,1).  This  principle  extends  to  any  class  of 
polynomials  orthogonal  on  the  real  interval  [a,b~\. 
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To  see  this,  choose  n  >  0  and  suppose  that  pn(x)  is  of  constant  sign  in  \_a,b~\.  Then 
<  p„(x) ,  p0(x)  >  t6  0,  which  contradicts  the  assumed  orthogonality.  Thus  by  continuity 
(and  the  Intermediate  Value  Theorem)  there  exists  a  zero  x}  e  (a,b). 

Suppose  that  .r,  is  a  double  root.  Then  p„(x)l(x  —  x,)2  would  be  a  polynomial  of  de- 
gree {n  —  2)  and  so 

0=  <pn(x),pn(x)l(x-xl)2>  =  <1  ,  (pn(x)l{x-x:))2>  >0 

which  is  a  contradiction.   Thus  the  zeros  are  simple. 

Now  suppose  that  pn(x)  has  exactly,/  zeros  xu  x2, ... ,  Xj  e  (a,b).  Then 

p„(x)(x  -  *,)(*  -  x2)  ...  (x  -  xj)  =  qn_j  (x)(x  -  xrfix  -  x2f  ...(x-  Xj)2 

where  qn_;  (x)  does  not  change  sign  in  (a,b)  and 

<  pn{x) ,  (x  -  xx){x  -  x2)  ...  (x  -  Xj)  >  =  <  qn_j  (x) ,  {x  -  xrfix  -  x2)2  ...{x-  Xj)2  >  . 

Since 

<  qn_j  (x)  ,  (x  -  xtfix  -  x2)2  ...(x-  Xj?  >  *  0 

and 

< Prix) ,  (x  -  xj(x  -  x2) ...  (x  -Xj)>  =0    forj  <  n, 

then  it  must  be  thaty>  n.  The  Fundamental  Theorem  of  Algebra  precludes^  n  and  so 
we  conclude^  n.  [Rcf.  5:  p.  236] 

Thus  all  the  zeros  of  p„(x)  are  real,  simple,  and  lie  in  the  interval  (a,b),  and  so  may 
be  ordered 

a  <  xl>n  <  x2n  <  -  <  xkn  <  -  <  xnn  <  b. 

The  interlacing  of  zeros  of  p„(x)  and  pn+l(x)  follows  from  the  Christolfel-Darboux  for- 
mula.   Recalling  that  the  leading  coefficient  kn  is  positive  for  all  pn(x),  then  (12)  gives 

( 1 3)  /?„'+]  (.y)  pn{x)  -  p'n{x)  pn+l{x)  >  0  ,    -oo  <  x  <  oo. 

Let  u  and  v  be  adjacent  zeros  of  pn(x).  Then 

(14)  p'n(u)p'n(v)<0 
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since  the  zeros  are  simple.  At  these  zeros,  inequality  (13)  reduces  to 

-/»/ViO<)>0 

and 

-p'n(v)Pn+i(v)>0- 
Multiply  these  two  inequalities  together  with  (14)  to  conclude 

pn+{(u)pn+l(v)<Q, 

so  /ViM  nas  a  zero  between  each  pair  of  consecutive  zeros  of  pn{x). 

Now  let  xn<n  denote  the  largest  zero  of  p„{x).  Observing  that  pn{x)  -*  oo  as  x  -»  oo, 
we  must  have  p'n{xn/)  >  0,  and  so  by  ( 13) 

Pn+l(xn/i)<°- 

But  pn+l{x)  -*  oo  as  x  ->  oo,  so  p„+l{x)  must  have  a  zero  to  the  right  of  x„^.  Similarly, 
pn+y{x)  must  have  a  zero  to  the  left  of  xiJX,  the  smallest  zero  ofpn(x).  Thus  all  n+l  zeros 
o[ pn+l(x)  are  accounted  for  and  interlace  those  of  p„(x)  [Ref.  8]: 

a  <  x\,n+\  <  x\,n  <  x2,n+\   <  '"  <  'Y^,/J+l  <  xk,n  <  xk+\,n+\  <  '"  <  xn,n+\  <  xn,n  <  xn+\,n+\  <  ^* 

E.     GENERATING  FUNCTIONS 

The  function  F(x,t)  having  a  formal  power  series  expansion  in  t 

oo 

«=o 

is  said  to  be  a  generating  function  for  the  set  \fn{x)}  [Ref.  9:  p.  129].   For  an  appropriately 
chosen  generating  function,  the  generated  set  {fn(x)}  is  a  class  of  orthogonal  polynomials. 
By  defining  a  class  in  this  way,  properties  of  the  polynomials  can  be  derived  from  the 
generating  function  itself.    . 
For  example,  consider 

oo 

(15)  F(x,t)  =  (1  -  2xt  +  /2)-|/2  =  £/,(*)/" 

>l=0 
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which  for  fixed  .y  is  the  Taylor  series  of  F(x,t)  centered  at  /  =  0.  If  we  restrict  x  to  the 
interval  [—1,1],  then  by  considering  the  singular  points  of  F(x,t),  we  may  conclude  that 
the  series  is  convergent  for  U|  <  1  [Rcf.  4:  p.  28J.  Thus  we  can  determine  f„(x)  for 
«  =  0,  1,2,...  by 

(16)  fn(x)  =  -Jr^r [(i  -  ixt  +  r2)",/2]:=o- 

3  1 

Equation  (16)  yields  f0(x)  =  1,    fl(x)=x,    f2(x)  =  —  x2  —  —  ,  etc.     We  also  note  from 

(16)  that 


=^[,!(i-.r-,],= 
=  i. 


We  now  derive  some  basic  properties  of  thcfn{x)  from  the  generating  function  in  (15). 
1.     Recurrence  Relation 

Differentiating  (15),  we  find 


oo 


(17)  -!f  =  ,(.-2.«+,J)-3,2  =  I/„'M<" 


oo 


(18)  -|^  =  (x  -  0(1  -  2xt  +  t2)-312  =  £  nfn{x)tn~\ 


n- 


Since  (x  —  t)  — t  —r—  =  0,  we  have 

ox  ot 


OO  '..  *  > 

(x  -  /)£/«  w  -  'Z  "/-W^1  -  ° 


n=l  /i=l 


which  becomes 


oo  oo  oo 


«=1  n—\ 
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Since  fa(x)  =  1,  we  can  start  the  sum  on  the  right  side  at  n  =  0,  then  re-index  so  that  it 
starts  at  n  =  1  again.   We  then  find 

oo  oo 

£  [  xfn'(x)  -  nfn(x)  ]  tn  =  £/„'_,  W  '"• 
«=]  n=l 

Equating  coefficients  off"  then  gives 

(19)  xfn\x)  -  nfn{x)  =/;_,  (x) ,    n>\. 

Rewrite  (17)  as 

oo 

(20)  (l-2.vr-f/2)-3/2  =  ^/„'(^"1,    t  +  0. 

Substituting  the  appropriate  expressions  from  (20),  (18),  and  (15),  respectively,  into  the 
identity 

(1  -  t2){l  -  2xt  +  t2)-312  ~  (2t)  (x  -  0(1  -  2xt  +  t2)~312  =  {l-2xt  +  tY12. 
now  gives 

oo  oo  oo 

(1  -  t^fn'W"-'  -  (2r)£  nfn(x)tn~l  =  £/n(*)r" 

n=\  «=1  n=0 

which,  upon  rearranging,  becomes 

oo  oo  oo  oo 

rt=l  rt=l  H=l  «=1 

Equating  coefficients  of /"and  gathering  terms  yields 

(2i)  (2n  +  !)/„(*)  =/;+1  (*)  -/„'_,  (*) ,  n  >  1. 

Substituting  (19)  into  (21)  gives 

which  by  a  shift  of  index  from  «  ->  «  —  1  becomes 

30 


(22)  xf-_,  {x)  =fn\x)  -  nfn_x{x) ,    n  >  2. 

Substituting  again  from  (19)  gives 

(23)  (*2  -  !)/„'(*)  =  nxfn(x)  -  nfn_{(x). 

Multiplying  (19)  through  by  (.v2  —  1),  we  have 

x{x2  -  \)fn'(x)  -  n(x2  -  \)fn{x)  =  (x2  -  !)/„'_,  (x). 

Now  substitute  for  (x2  —  l)fn'(x)  and  {x2  —  1)/,!,  (x)  from  (23)  and  gather  terms  to  get 

(24)  nfn(x)  =  {In  -  1)  */„_,(*)  -  (n  -  l)fn_2(x) ,    n>2. 

Equation  (24)  is  a  three-term  recurrence  relation  for  [fn(x)}.  The  advantage  of 
this  form  is  that  beginning  with  f^x)  and  f^x),  we  can  now  generate  any  member  of 
{fn(x)}  by  iterating  (24)  and  thus  avoid  the  differentiation  in  (16).  Note  that  since 
fo(x)  =  1  -  fi(x)  =  *.  (24)  implies  that  [fn{x)}  is  a  set  of  polynomials.   [Ref.  9:  pp.  159-160] 

2.  Ordinary  Differential  Equation 

We  continue  the  same  line  of  reasoning  to  extract  additional  information  about 
[fn{x)}.  Differentiating  (22)  yields 

(25)  x/;_x  (x)  =f„\x)  -  in  +  1  )/„•_,  (*). 

From  (19)  we  have  /,'_,  (x)  and,  after  differentiating,  f"n_x  (x).   Substituting  these  ex- 
pressions into  (25)  gives 

x  [xfn'(x)  -(n-  \)fn'(x)]  =/„"(*)  -(n+  l)[_xfn'(x)  -  «/„(*)] 
which  when  rearranged  becomes  a  second  order  ordinary  differential  equation 

(26)  (1  -  x2)fn'(x)  -  2xfn'(x)  +  n(n  +  l)fn(x)  =  0. 

The  lfn(x)}  are  solutions  of  (26)  for  n  =  0,  1,  2, ... .  [Ref.  9:  pp.  160- 161] 

3.  Orthogonality 

Rewriting  (26)  as 

(27)  -£■  [(1  -  x2)fn'(x)]  +  n(n  +  !)/„(*)  =  0 
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we  recognize  the  structure  of  a  Sturm-Liouville  eigenvalue  problem.    Since  the  points 
„r  =  +  1  are  singular  points,  we  require  ihatfn(x)  andfn'(x)  be  finite  asx-*+  1.  From  the 
associated  theory  of  the  singular  Sturm-Liouville  problem,  we  conclude  that  the  {fjx)} 
are  orthogonal  on  the  interval  [  — 1,1]  with  weight  function  w(x)  —  1. 
Alternately,  we  combine  (26)  with 


dx 
where  n  ¥=  m,  to  obtain 


[(1  -  **)/„'(*)]  +  m{m  +  \)fjx)  =  0 


(28)  UX)  ~d7  ^  ~  x2)f"'W]  "/-(*)  ~h  [('  "  X2)frn'{*)} 

+  0("  +  1)  -  m(m  +  l)]fn(x)fjx)  =  0. 

Since 


-£-  [(1  -  x2)lfm(x)fn'(x)  -fm'(x)fn(x)}] 

=  -^[^W{(l-^)/n'W}  -fn(x){(l-x2)fm'(x)}] 
-/«(*)  £  [0  -  *2)//W]  "/«W  ■£"  [(1  "  *2)/m'W], 

we  can  write  (28)  as 

[(1  -  x2){fm(x)fn'(x)  -fm'(x)fn(x)}]  +  [n2  -m2  +  n-  m]fn(x)fm(x)  =  0 


dx 


or 


(n  -  m){n  +  m  +  \)fn{x)fm{x)  =  -£-  [(l  -  *2){/mW„(*)  -/„,(*)/„'(*)}]• 
Integrating  from  x  =  —  1  to  x=  1,  we  obtain 

(„  _  m)(„  +  w  +  DJ1  fn{x)fm{x)dx  =  [(1  -  *2){/mW„to  "/«(*)/«'(*)}]_!, 
Since  (1  —  .v2)  =  0  at  jc  =  ±  1,  we  have 
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(„  _  m)(„  +  m  +  n  J    fn{x)fjx)dx  =  0. 
Recalling  that  «  =£  m,  we  conclude 

[  Mx)fm{x)dx  =  0 

J-i 

i.e.,  the  polynomials  (/i(jc)}  are  orthogonal  on  the  interval  [—1,1]  with  respect  to  the 
weight  function  w(x)  =  1.  [Ref.  9:  pp.  173-174] 

These  are  the  Legendre  polynomials  that  were  previously  defined  using  the 
Gram-Schmidt  process.   Thus  ihef„(x)  of  (15)  are  in  fact  Pn(x)  and  (15)  may  be  written 

oo 

(29)  (i-2xt  +  t2)-]l2  =  ^Pn(x)tn, 

establishing  the  equivalence  of  the  generating  function  definition  with  the  Gram-Schmidt 
definition. 

Legendre  and  Laplace  concluded  that  the  Pn(x)  in  (29)  were  polynomials  of  de- 
gree n   in  the  variable  x  by  examining  a  series  expansion  of  the  function 

(l-lxr+r2)-111  -  7  rn  7     (~°,  ^T    {2x)n-2m. 
v  I  Lj       Lu        m\  (n-2m)\ 

They  reasoned  the  orthogonality  directly.    From  (29)  we  write 

oo 

(l-2jcr+r2)~1/2  -  Yj  Pfa)  r 


i=0 


and 


(l-2xs+52)~,/2  =  ]T /y.v)/ 


y=0 

Multiplying  these  power  series  together  via  the  Cauchy  product  formula  yields 
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yj  1—  2xr+r    yj l—2xs+s* 


= Z  Z **w  p<^w rk  r 


«=o  /<=o 

oo 


«=0 


=  Z    Z  pM  pM) 


i+j=n 


r'/ 


Integrating  from  -1  to  1  with  respect  to  x  gives 


I 


dx 


/  2        /  2 

1    J  1  —  2.rr  +  r     w  1  —  2xs  +  5         «=o  j+y=« 


oo  I —      . 

=  Z  Z    J    pM  pM) dx 


■V. 


Through  tedious  calculation,  the  left-hand  side  of  this  expression  becomes 


1       ,       1  + V™ 
log 


Jrs  '  1  -  JTs 


a  function  of  the  product  rs.     From  this  we  conclude  that  the  coefficients  of  the  terms 
in  the  series  on  the  right-hand  side  are  zero  when  /  andy   differ,  i.e., 


r  PMPjddx-o,  i+j 

J-l 


and  the  orthogonality  is  established.    Finally, 


(\-2xr+r) 


2\  — 1/2 


oo 

-TV 


*-1        \-r       Z-, 


or  P„(l)  =  1  and  so  the  Pn(x)  are  in  fact  Legendre  polynomials.   [Ref.  8] 

In  a  similar  fashion,  the  norm  of  the  Legendre  polynomials  can  also  be  obtained 
from  the  generating  function.    Let 


;„  =  /!;,=J'  [Pn(x)fdx. 


Square  both  sides  of  (29)  and  integrate  with  respect  to  x   on  [  —1,1] 
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oo       oo       /        .  \ 

E  E  (  J      Pm(x)  Pn(x)  dx     tm+n  =  J 

«J=0  /!=0     \       _1  / 


* T  dx. 

-i     1  -  2*/  +  r 

By  orthogonality,  the  left-hand  side  is  zero  unless  m  =  n,  while  the  right-hand  side  is 
integrable  in  closed  form: 


-l 


/      \  l-/ 

This  function  can  be  expressed  as  a  difference  of  two  logarithms,  each  of  which  has  a 
convergent  Maclaurin  series  expansion  in  (  —1,1).  When  combined  this  yields 


S^-Z  ■= 


+  1 

n=0  n=0 


t2n. 


2 
Comparing  coefficients  of/2"  on  both  sides  gives  c„  =  — —  . 

Since  there  is  no  systematic  theory  for  determining  generating  functions,  finding 
one  for  a  polynomial  class  can  be  a  problem.  The  work  above  bears  this  out.  Unfor- 
tunately, the  proofs  above  do  not  easily  generalize  to  related  classes.  With  this  in  mind, 
let  us  summarize  the  key  steps  in  proving  that  (15)  is  a  generating  function  for  the 
Legendre  polynomials.  First  we  established  that^(l)=  1  for  n  >  0.  Next,  we  showed 
that 


n=0 


for  fn(x)  a  polynomial  of  degree  n    in  the  variable  x.  Finally,  we  showed  that  these 
polynomials  were  orthogonal  on  the  interval  [—1,1] 


(x)  dx  =  0  ,  m^  n. 


These  three  points  are  sufficient  to  show  that  the  generated  class  is  the  Legendre  class 
of  orthogonal  polynomials. 
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We  now  provide  an  alternate  proof  attributed  to  Ilcrmitc.  Our  interest  in  this 
proof  is  mainly  the  technique  which  suggests  a  method  of  generalization  that  we  will 
take  advantage  of  in  the  next  chapter. 

We  begin  with  (15)  reproduced  below 

oo 

Multiply  both  sides  of  this  equation  by  x*and  integrate  from  -1  to  1  to  obtain 


(30)  I"'  x<dx  -£  tn  [  xkfn(x)  dx. 

-1    Jl-Zxt+t2        n=0        -1 


Now  change  variables  from  x  to  y  via 

(1-2^  +  ,2)-,/2=1-/k 

giving 

t(l-y2) 
x  = +  y  ,    dx  =  ( 1  -  ty)dy. 

The  left-hand  side  of  (30)  becomes 

I,  (i-o»  {i-'y)dy 

or  simply 

(;[i<(.-/)+>]"4. 

Expanding  the  integrand  of  this  last  expression  via  the  Binomial  Theorem,  we  obtain 


y-o 
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When  written  as 


we  can  identify  the  form  as  a  polynomial  of  degree  k  in  the  variable  t.    Comparing  this 
form  to  the  right-hand  side  of  (30),  we  conclude  that 


x  fn(x)  dx  =  0  ,    for  n  >  k, 
J-1 


i.  Q.,f„{x)  is  orthogonal  to  all  polynomials  of  lower  degree.  The  same  argument  as  before 
gives  7^(1)  =  1  for  n  >  0  and  the  proof  is  complete. 

F.     HYPERGEOMETRIC  SERIES 

The  term  "hypergeometric"  was  used  in  1655  to  distinguish  a  series  that  was  "be- 
yond" the  ordinary  geometric  series  1  +  x  +  x2  +  ••• .  In  1812,  Gauss  presented  the  power 
scries 

ab_^_      a(a+l)b(b+l)   x2       a(a+  \)(a+2)b(b+  \)(b  +  2)    x2 
+    c     V.  +         c(c+l)  2!   +  c(c+l)(c  +  2)  3!   + '" 

c  =£  0,  —  1,  —2,  ...  which  is  known  as  Gauss'  series  or  the    ordinary  hypergeometric  series 
[Ref.  10]. 

Convergence  of  this  series  for  |„r|  <  1  follows  directly  from  the  Ratio  test.  By 
Raabe's  test,  convergence  can  be  shown  for  |.v|  =1  when  (c  —  a  —  b)  >  0  [Ref.  11:  p.  5]. 
Gauss  also  introduced  the  notation  2F{[_a,b ;  c ;  .v]  for  this  series.  Note  that 
1Fl\_a,b  ;  c  ;  x~\  may  be  considered  as  much  a  function  of  four  variables  as  a  series  in  x. 
[Ref.  12:  p.  1] 

With  the  shifted  factorial,  the  ordinary  hypergeometric  series  can  be  expressed 


tt=0 
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Below  are  some  examples  of  important  functions  which  can  be  expressed  as  ordinary 
hypergeomctric  series. 

Example  1:  log(l  +  x)  =  .x-2/7,[l,  1  ;  2  ;  —  x] 

Example  2:  siir'(-v)  =  x2F1l\/2,  1/2  ;  3/2  ;  *2] 

Example  3:  tarr'(jc)  =  jc  ^[1/2,  1  ;  3/2  ;  -x2] 

The  generalized  hypergeometric  series  is  formed  by  extending  the  number  of  pa- 
rameters, an  idea  attributed  to  Clausen  [Ref.  1 1:  p.  40]. 


oo 

ZMnMn  ■■■(ar)n         n 
{b,)n{b2)n...(bs)nn\  X  ■ 


rFs\_ax,  a2,  ...  ,ar\blt  b2,  ...  ,  bs ;  .r]  = 

Note  that  since  {a)n^j{a)n  =  n  +  a,  a  hypergeometric  series  £  cn  •*"  is  characterized  by  the 
fact  that  the  ratio  cn+Jcn  of  coefficients  is  a  quotient  of  two  polynomials  in  the  index  /;, 
i.e.,  a  rational  function  of  n. 

The  Ratio  test  can  be  used  to  show  convergence  for  all  values  of  x  when  r  <  s  and 
for  |.x|  <  1  when  r  =  s+  1.  When  r  >  s  +  1,  the  series  diverges  for  all  x  i=-  0  and  the 
function  is  defined  only  if  the  series  terminates.  The  series  terminates  when  one  or  more 
of  the  numerator  parameters  a,  is  zero  or  a  negative  integer  [Ref.  11:  p.  45].  This  is  an 
important  characteristic  of  the  hypergeometric  series  that  will  be  used  later.  A  power 
series  that  terminates  gives  a  polynomial  which  is  defined  for  all  x.  In  this  case,  the  pa- 
rameters bu  ...  ,  bs  may  be  negative  integers  as  long  as  the  series  terminates  before  a  zero 
is  introduced  into  a  denominator  term. 

Examples  of  the  generalized  hypergeometric  series  include  familiar  functions  such 


as: 


Example  4:  (1  +  x)°  =  ,F0[  -a ;  —  ;  -jc] 
Example  5:  e'  =  0F0[  —  ;  —  ;  x] 
Example  6:  sin  x  =  x  0F{ [  —  ;  3/2  ;  —  x2/4] 
Example  7:  cos  x  =  0F,[  —  ;  1/2  ;  —  *2/4] 
Example  8:   The  Bessel  function  of  order  a 

.       (*/2>Vi[-  ;«  +  l;  -*2/4] 

J<{x)  = nTTT) 

where  dashes  indicate  the  absence  of  parameters,  i.e.,  when  r  =  0  or  5  =  0.  We  will  also 
use  a  common  alternate  notation 
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rfs 


flj,  ...  ,  ar 
_&,,...  ,bs 


,x 


for  either  the  ordinary  or  generalized  hypergeometric  series.   [Rcf.  12:  p.  4] 
1.     Chu-Vandermonde  Sum 
The  Chu-Vandermonde  sum 


2F\l-n,  a;c;  1]  = 


(c  -  a)» 
(c)n 


is  one  of  many  useful  summation  formulas.  Since  this  one  will  be  used  in  a  later  chapter, 
the  proof  is  provided  below. 

Basically,  this  is  a  consequence  of  the  General  Binomial  theorem 


(1-*)' 


oo 


)k      Ic 


k=0 


Starting  with  the  identity  (1  —  x)~a(\  —  x)-6  =  (1  —  x)-'-b,  expand  both  sides.  Us- 
ing the  Cauchy  product  on  the  left  side  and  the  General  Binomial  theorem  on  the  right, 
we  have 


n=0 


n=0 


(f  +  b)n     „ 
0)*      *' 


"     {a)k{b)„-k     -  . 

where  c„  =  X    ,.^  ,  ^ —  .  Equating  coefiicients  ol  xn,  we  have 


*=o    (1)*(1)„- 


z 

k=0 


(a)k(b)n_k        (a  +  b)„ 


(1  )*(!)«-*  (1)« 


In  order  to  express  the  left  side  as  an  ordinary  hypergeometric  series,  2Fl,  mul- 
tiply both  sides  by  (1)„  and  use  the  identity  (I  )„/(!)„_*  =  (  — 1)*(  —  n)k  to  obtain 


(31) 


n 

z 

k=0 


( —n)k(a)h  u 

(-lf(b)n_k  =  (a  +  h)n. 


(l)i 


Next, 
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(-iA^_,=(-D"(-ir^_, 

=  (-l)>\-lf-k(b)(b+[)...(b  +  n-k-l) 
=  {-\)"{-b){-b-  l)...(-b-n  +  k+l) 
n{-b-n  +  1)„ 


(-b-n+l),, 
(-b-n+l). 


=  (-D 


(-b~n+l)k 


Using  the  above  result,  equation  (31)  becomes 


n 


( ~n)M)k 


^         (1),         {-b-n+\)k 


=  (-D' 


(a  +  b)n 


{-b-n+l)n 


Let  c  =  —b  —  n  +  \,  and  substitute  to  get 


2Fil-n,  a;c;  1]  =  (-1) 


„  (-c  +  a-n+  \)n 


(c)n 


(-!)"( -c  + a- w  +  l)...(-c  +  «) 

(C)n 

(c  —  a  +  n  —  1)  ...  (c  —  a) 


{C)n 


{C)n 
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IV.     JACOB!  POLYNOMIALS  AND  SPECIAL  CASES 

This  chapter  focuses  on  the  orthogonal  class  known  as  Jacobi  polynomials,  P^'p\x). 
Until  recently,  the  classes  of  orthogonal  polynomials  considered  "classical"  were  usually 
given  to  be  Jacobi,  Gegenbauer  (also  called  ultraspherical),  Chebyshev  (of  first  and  sec- 
ond kind),  Legendre  (also  called  spherical),  Laguerre,  and  Hermite.  The  Jacobi 
polynomials  hold  a  key  position  in  this  list  since  the  remaining  classes  can  be  viewed  as 
special  or  limiting  cases  of  this  class.  Today,  the  classical  orthogonal  polynomials  are 
taken  to  be  special  or  limiting  cases  of  either  of  two  very  general  orthogonal  classes 
known  as  the  Askey-W  ilson  polynomials  and  the  q-Racah  polynomials,  between  which 
we  will  establish  a  formal  equivalence.  Because  of  their  complexity,  these  classes  are 
described  in  Chapter  VI  after  the  necessary  additional  theory  has  been  developed. 

The  results  derived  in  the  text  that  follows  are  arranged  in  tabular  form  by  class  at 
the  end  of  this  chapter. 

A.     JACOBI  POLYNOMIALS 

1.     Definition  /  Orthogonality 

The  Jacobi  polynomials  P(*'P)(x)  are  generated  by  applying  the  orthogonalization 
step  of  the  Gram-Schmidt  process  to  the  standard  basis  {l,x, x2,  ...  }  of  P[  —1,1],  with 
respect  to  the  weight  function  given  by  a  continuous  beta  distribution  on  [  — 1,1] 


w(x;oc,  /?)  =  (!  -xf{l  +  x) 


i< 


for  a  >  —  1,  /?  >  —  1,  i.e., 


| '   Pt  %)  Fj*  "'(.V)  (1  -  xf  (1  +  xf  dx  =  [/■<«'  ">]- 'im.n- 


The  Jacobi  polynomials  can  also  be  represented  by  hypergcometric  series 


(i) 


^.«w=  («+!)„ 


n\ 


r     2^1 


— nf  n  +  a  +  /?  +  1     [  —  x 
a+1  ;      2 


This  set  of  polynomials  is  thus  standardized  (as  was  done  for  the  Legendre  class  in 
Chapter  II,  Section  B)  : 
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e'*(l)= 


n  +  oc\       (a+  1)„ 


A/! 


We  shall  demonstrate  that  these  two  characterizations  of  the  Jacobi  polynomials  are 
indeed  equivalent. 

To  show  that  the  polynomials  defined  in  (1)  are  in  fact  orthogonal  with  respect 
to  w(x  ;  a,  /?)  =  (1  -  jc)"(1  +  xf  on  [  — 1,1],  it  suffices  to  show  that  P^'p\x)  is  orthogonal 
to  one  polynomial  of  degree  m  for  m  =  0,  1, 2, ...,«—  1.  This  is  because  any  /'^(.v), 
0  <;  m  <  n  —  1,  can  be  expressed  as  a  linear  combination  of  such  polynomials.  While  any 
polynomial  of  degree  m  could  be  used,  we  choose  (1  +  x)m  for  reasons  that  will  become 
apparent. 

To  establish  orthogonality,  we  consider 


P^\x),(l+x)m> 


P<a'  %)  (1+^(1-^(1+  xf  dx. 


-1 


By  (1),  the  right-hand  side  becomes 


(-n)k(n  +  a  +  P+\)k    fi 


k=0 


(«+  i)kk\r 


J" 


\w+/? 


Using  the  last  result  of  Chapter  II,  Section  C.l,  the  change  of  variable  x=  I  —2t  yields 


(«  +  0„ 


n\ 


Jt=0 


(a+  1)*A!2* 


(a  +  1)„ 


«! 


/<=0 


(-w)fe(w +  «  +  /?  +  !)„  \    w+flt^+1       rXM-a+  l)r(/»  +  /?+  1) 

A!  J  (a+l)ir(Hw  +  a  +  /?  +  2) 


Identities  for  the  gamma  function  allow  us  to  simplify  this  expression  to  the  form 
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23"f/7+m+1  r(a  +7i  +  1)  Vim  +  ft  +  1)  V  (  ~n)k  ("  +  «  +  /?  +  1)/, 


T(m  +  a  +  /i  +  2) 


(m  +  a  +  /?  +  2)fc  /<! 


/c=0 


which  can  then  be  written 


2a+/?+"7+1r(g  +  /7+l)r(m  +  /?+)) 
f(m  +  a  +  /J  +  2) 


2*i 


■«  ,  w  +  a  +  ft  +  1 
m  +  a  +  /?  +  2 


;i 


The  Chu-Vandermonde  sum  from  Chapter  II  allows  us  to  write  this  as 
2a+p+m+\  r(a  +  n+{)  T{rn  +  0+1)         (m  +  1  -  n)n 


T{m  +  a  +  ft  +  2) 


{m  +  a  +  ft  +  2)n 


Thus, 


<P^%),(l+xf>  = 


0,  m  =  0,  1,2, ...  ,n-  1 


where 

n\  2n+*+f]+l  r(n  +  a  +  1)  Yin  +  ft  +  1) 
C"  ~  r(2/i  +  a  +  ft  +  2) 

which  justifies  the  orthogonality.  It  is  possible  to  extract  the  value  of  [/£'">]~1  by  modi- 
fying this  argument,  but  we  defer  this  computation  till  the  next  section,  where  it  will  be 
easier. 

2.     Ordinary  Differential  Equation  /  Rodrigues'  Formula  /  Norm 

We  begin  deriving  the  ordinary  differential  equation  for  the  Jacobi  polynomials 
by  noting  the  general  formal  result: 


(2) 


d      p 
dx 


lA»-  ,b 


\x 

s       J 


a,  ...  nr 
*.  -  bs 


a,  +  1,  ...  ,  ar  +  1 
_/>!  +  1,...  ,bs+  I 


;-v 


This  can  be  seen  from  the  definition  of  the  generalized  hypergeometric  series  by 


writing 
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d      p 
dx  r  s 


«,, ...  ,ar 
L/>,,...  ,b 


,x 

s        J 


d 

dx 


oo 


I 


(a])k...(ar)k      xk 


k=0 

(a])k...(ar)k 


£-<    {bx)k...{bs)k      [k-\)\ 


-z 

fe=0 


Noting  that  (a)*+1  =  a  (a  +  1)*,  we  have 


dx   r  s 


au  ... ,  ar 
bu...,b: 


Lu\ 


S  J 


a,  ...  a, 


6,  ...6 


ou 


(a,  +  1)^.  ..{ar+  l)k 


x 
(bl  +  l)k...(bs+\)k      k\ 


from  which  (2)  follows.    Differentiation  of  the  hypergeometric  scries  is  justified  by  re- 
calling that  the  power  series  is  convergent  for  |  jc  |  <  1  when  r  =  s  +  1. 
We  apply  this  result  to  the  Jacobi  polynomials  to  obtain 


<L  *(*•*><*        d    /^(«+D. 


dx   ?n 


V) 


dx 


-Z7, 


//! 


— n,  n  +  a.  +  ft  +  1      l  —  x 

a  +  1  :  ~T~ 


_1   \    (a+l)„      -n(n  +  a  +  fl+\) 


2^1 


2    /         «!  a+1 

(tt  +  a  +  /?+l)(a  +  2)„_,        («-!)!         (a+li/?+1) 


-/»+lfH  +  a  +  /?  +  2 

a  +  2 


—  jc 


2(»-l)! 


(«  +  2)„_i 


n-\ 


w 


which  simplifies  to 


(3) 


*  i 


rfx 


W  = 


(AH-a+_gj-l)      (a+1|  ^+i; 


n-\ 


(X). 


With  these  results  established,  we  now  consider  the  orthogonality  property  of 
the  class:    set 


(4) 


J-l 


By  (3),  we  can  write 
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Integration  by  parts  yields 

noting  that  the  boundary  terms  vanish,  so  that 

(6)  ^n    =   ~  +~2+  g    \  \  Pfcl'  *-%)  1n+l  W  (1  "  *)- '  (1  +  */"'  * 

where 

(7)  ^+1W  =  (1  -  *2)  -£■  p^(x)  -  a  (1  +  x)  P(;-%)  +  fi{l~x)  Pn'%) 

is  a  polynomial  of  degree  (n  +  1).  Thus  by  (6)  we  may  express  qn+[(x)  as  a  linear  combi- 
nation of  Jacobi  polynomials 

qn^x)  =  Y,CjPf '?-%). 

7=0 

We  would  now  like  to  show  that  c0  =  c{  =  ■■■  =  c„  =  0,  so  that  only  the  last  term 
of  the  summation  survives.  The  constants  c;,  j  =  0,  1,2, ...  , n  can  be  determined  by 
substituting  this  expression  for  qn+l(x)  into  (6)  and  using  (4). 


#f  =     Pt  %)  i%  %)  ( i  -  *T  ( i  +  *f  dx 

ttj SCJ'  f«'Hlwf'"l'H,w(i-'r,(i+/-1*. 


m  + 

7=0 


For  each  m<n—  1,  by  (4)  the  left-hand  side  /m„  =  0,  but  the  right-hand  side  is 
zero  unless  j  =  m  +  1  <  n.  Thus 
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Cj  =  0  ,  j  <  n 


which  means  from  (6) 


and 

(9)  ^+1W  =  c,+lJpt"/'/?"1)W. 

To  determine  the  remaining  constant  cn+u  we  note  that  by  (7)  and  (9) 

(l-*2)  -£  Pf%)  -  a  (1+,)  Pt%)  +  *(!-*)  if '«(*)  =  cn+1  /»$■  »-%). 

Letting  jc=  1,  then  from  the  hypergeometric  series  definition  of  the  Jacobi  polynomials 
we  have  I**' p\l)  =  —    ,       which  gives 

n\  C"+]    (n+i)\ 

and  so 

(10)  C|I+1--2(/f+l). 

Combining  (5),  (6),  (9),  and  (10),  we  have 

(id  ■£■[ p<?>%) d-xTd+xf  ]  =  -2 (n+\) [ />£-/■  '-»(,) d-.r1  d+xf-1]. 

To  obtain  the  second  order  ordinary  differential  equation,  change  n-*n—l,  a->a+l, 
and  P  ->/?+!  in  (11)  to  give 


ah 

Then  by  (3), 


f  [/ft1'  ^'VKl-xf1  (l+^+1  ]  =  -In  Pj*»(x)  (l-xfd+xf. 
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d 

dx 


n+tx+(]+\       dx 


Lp^\x)d-xri(i^f+] 


><«,  p) 


=  -2nPr,"(x)(l-xr(l+xy, 


i.  e.,>'=  P*   \x)  satisfies 


d 
dx 


('^c^'i 


=  -n  (n+a+fl+l)  (1-jc)"  {\+xfy. 


Using  the  product  rule  to  expand  the  left-hand  side,  we  have 

(l-xTd+xf  [-(a+1)  (1+Jf)y  +  GJ+1)  (l-x)y'  +  (l-Jf2)/  ] 
=  -/i  («+«+/?+ 1)  (l-jcf  (1+jc)% 

which,  when  simplified,  becomes  the  second  order  ordinary  differential  equation  for  the 
Jacobi  polynomials  j>  =  P(*'P)(x)  : 


(12) 


(l-jc2)/'  +  [  ifi-a)  -  (a+P+2)  x]y  +  n  (n+z+p+l)y  =  0. 


The  reader  is  invited  to  compare  this  result  with  the  second  order  ordinary  differential 
equation  for  the  Chebyshev  polynomials  {Tn(x)}  given  by  Equation  (5)  in  Chapter  I. 
By  iterating  (11)/:   times,  we  obtain 


<t 


dx 


[  pW\x)  (1-xf  (l+xf  ]  =  (-l)k  2k  (n+l)k  [  />&*  P-k\x)(l-xrk(l+xf-k  ]• 


Setting  a  -*  v.+k,  and  /?  -*  fi+k  gives 


(l-xT{l+xfpfrk%)  = 


(-1) 


2"  («+!)*      dx 


i-  [  ( ,  _*)«+*  ( i  +xf+*  pfa*.  m\x)  ]f 


or  equivalently  in  terms  of  the  weight  function, 


>(«.  /?V 


(-1 


<T 


w(,;«,  /*)/>£/'(*)  = 

2   (w+lj/j      tf.x- 


[  »K*  J  «+M+*)  #*  WM  ], 


the  general  Rodrigues  formula    for  the  Jacobi  polynomials.    Letting  n  =  0  and  noting 
P(;'p\x)=  1  gives 


A  p(*.  /?), 


(i-.vr(i+xy>/^',,u)  =  — 


(-If    / 


2"  A!       «rfr* 


[(!-,)•+*  (1+jc)^*], 
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or 

w(x  ;  a,  (I)  /£'  %)  -  -tli-   -£-  [  w(*  ;  a+*,  0+A)  ], 

2A  /<!       <£c* 

the  classical  Rodrigues  formula  for  the  Jacobi  polynomials. 

Finally,  we  can  use  these  ideas  to  obtain  the  value  of  [h(*'P)Yl.  By  direct  com- 
putation, 


J-i 


2a+/?+1r(a+l)r(/?4-l) 
r(a+^+2) 

We  will  use  this  below.   Combining  (8),  (10),  and  (4)  when  m  =  n  yields 

[,(*, /?)!-!         4(w+l)     r    (,-1,/J-in-i 

L^     J    -  „+a+/J  LAn-i       J   • 

Making  the  changes  n  -*  a?— 1,  a  -»  a+1,  and  /?  -*  /?+!,  this  may  be  rewritten  as 

KtT  =  "T'K''ft"]"' 

Iterating  this  relation  n   times  produces 

r(2w+a+/?+l)         22"+a+/?+1  r(/i+a+l)  !>+/?+ 1) 


22%!  r(/7+a+/?+l)  r(2«+a+/?+2) 

2a+/?+1  r(At+a+l)r(A;+/?+l) 
"  (2«+a+/?+l)  «!  r(«+a+/?+l)  * 

3.     Generating  Function 

As  was  mentioned  in  Chapter  III,  finding  a  generating  function  for  a  class  of 
orthogonal  polynomials  can  be  a  challenging  task.  Fortunately,  a  generating  function 
for  the  Jacobi  polynomials  can  be  found  by  mimicking  a  technique  attributed  to  Hcrmite 
in  his  work  with  the  Legcndrc  polynomials.  (Sec  Chapter  III,  Section  F.3.) 

We  seek  a  generating  function 
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(13) 


rJT> 


n=Q 


for  constants  cn ,  n  =  0,  1,  2, ... ,  with 


j 


,(«,/?),  ,    k 


P^'p>{x)xK(l-xT(l+xfdx  =  0,  k<n. 


P  ^^ 


To  begin,  we  multiply  (13)  by  xk  (1  —  x)*  (1  -f  ^  and  integrate  from  -1  to  1  to  obtain 


(14)       .v/<JF(r,x)(l-x)a(l+^^=yc„ 
J-1  »=n 


;?=0       L 


pr  V)*  (l-.vfo  +*f  </„* 


-1 


Note  that  the  summation  of  the  right-hand  side  is  from  0  to  k  because  the  orthogonality 
property  makes  each  term  zero  for  n  >  k.  We  observe  that  the  unknowns  to  be  deter- 
mined in  (14)  are  F{r,x)  and  cn.  When  one  is  given,  the  other  can  be  found,  so  we  first 
consider  the  left-hand  side  of  (14) 


[    xkF{r.x)(\  -.vf  (1  +  xfdx. 


i /  i-y2\ 

Setting  y/l—  2xr+r2  =  \—ry  and  substituting  for  x  =  y  +  I  — - —  I  r  yields 


-l  L 


y+ 


(lV) 


F  [r,  y+ 


(1-/) 


l-y- 


(lV) 


l+y+ 


(lV) 


{l-ry)dy. 


Factoring  out  the  terms  (1— y)*  and  {\+y)p,  we  have 


(15) 


I 


v+ 


OV) 


(1-/) 

[r<y+— -y- 

\+y)r  ' 


1- 


1  + 


(\-v)r 


\—ry)dy. 


We  would  now  like  to  make  a  judicious  choice  of/7  in  order  to  facilitate  calcu- 
lations. Following  the  lead  of  the  Legcndre  polynomials  (Chapter  III.  Section  F.3) 
suppose  F  is  such  that 


49 


(16)  F[r,y+ 


OV) 


{\+y)r 


1+ 


(1-v 


(\-ry)=l, 


then  the  integral  in  (15)  becomes 


f> 


(i-/) 


2\        -, 


(l-y)a  (l+yf  dy. 


We  note  that  the  integrand  is  a  polynomial  of  degree  k    in  the  variable  r   as  is 
the  right-hand  side  of  (14).    From  (16),  we  get 


F  [r,y+ 


(1-/) 


=  {\-ry) 


-l 


d+y)r 


1+ 


d-j') 


which    becomes,    via    1—  2xr+r2  =  (1—  ry)2,    the    generating    function    for    the    Jacobi 
polynomials: 

F(r,x)  =  2a+/?  (l-2xr+r2)-1'2  [l-,+  (l-2.rr+,2)1/2]-a  [\+r+(l-2xr+r2y12]^. 
Using  this  generating  function  in  (13)  gives 


£  cn  P^'%)  rn  =  2a+/?  (l-2.rr+,2)-1/2  [l-r+(l-2^+,2),/2]-*  [l+r+(l-2^+r2),/2]-/? 

from  which  we  can  determine  cn  by  setting  x=  1  and  recalling  that  P(*,p\\)  = r—2-. 

Thus, 


Z  c»  "^T^  '"  =  2^  (l-r)-'  [2"*  (1-rpJ  2"/?  =  (l-rp"1 


which  bv  the  Binomial  Theorem  becomes 


oo 

E 

«=0 


c« i r    =  /  .  i r 

n=0 


and  so  c„  =  1  for  all  /?.  [Ref.  8] 
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B.     SPECIAL  AND  LIMITING  CASES 

With  the  structure  of  the  Jacob!  polynomials  established,  we  turn  now  to  the  role 
of  the  parameters  a  and  /?. 
1.     Special  Cases 

For  certain  choices  of  the  parameters  a  and  /?,  we  find  that  the  classes  previously 
examined  and  several  new  classes  are  produced  as  special  cases  of  the  Jacobi  class. 
These  subclasses  inherit  the  structure  of  the  parent  class  which  often  provides  a  direct 
way  to  establish  specific  properties  (i.e.,  polynomial  nature,  orthogonality,  etc.) 
Table  1  provides  the  choice  of  parameters  for  selected  classes. 

Table   1.     PARAMETERS  FOR  SELECTED  CLASSES 


Class 

Parameters 

Jacobi 

<x>-l,    /?>-l 

Gegenbauer 

a  =  (]  =  A-\j2 

Chebyshev, 
First  Kind 

ct  =  /?  =  -l/2 

Chebyshev, 
Second  Kind 

«  =  /?  =  1/2 

Legendre 

a  =  fi  =  0 

2.     Limiting  Cases 

In  this  section,  we  briefly  examine  the  Laguerre  polynomials  and  the  Ilermite 
polynomials.    Using  the  hypergeometric  series  definition  of  the  Jacobi  polynomials,  we 
show  how  the  Laguerre  class  is  a  limiting  case  of  the  Jacobi  class. 
a.     Laguerre  Polynomials 

The  Laguerre  polynomials,  V*\x)t  are  defined  in  terms  of  hypergeometric 
series 


(17) 


M 


TO  = 


n+cc 


—n 


a+1 


\x 


The  first  relationship  to  establish  is 
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tfW-Jmjfr-fl-f 


To  show  this,  we  begin  with  the  hypergeometric  series  representation  for  the  Jacobi 
polynomials 


I  "+a  l  .F, 


/*••  «>W  = 


A/ 


-«,  «+<*+/?+ 1       1-v 
a+1  ;  ~2~ 


2r 
When  y  =  1 t-  .  we  nave 


n+a 


iFi 


■n,  n+a+p+l      x 
a+1  ;  7 


Writing  out  the  power  series,  we  obtain 


A7  +  0C 

«+a 

A7 


E 

oo 

z 

A:=0 


(-/f^  (/!+«+/?+!)*      ^ 

(a+1)*  A!  /** 

(-«)***      («+a+/?+l)A 


(a+1)*/:!  ^ 


In  the  limit  as  /?  -+  oo,  the  ratio 


(*+«+/?+ l)fe 


and  so 


completing  the  proof.   [Ref.  13:  p.  103] 

To  obtain  the  orthogonality  relation  for  the  Laguerre  polynomials,  we  start 
with  the  orthogonality  relation  for  the  Jacobi  polynomials 

f '  ft  P)i  y)  Pt  /?)( y)  ( i  -yf  ( I  +yf  dy  =  o,    m  *  «. 

J-i 
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Letting  y  =  1 — -r-  gives 


P£  »( l-JZ. )  P<«.  » 


1- 


P      \  P 


>.x    \* 


"-?) 


/?     2 


^/x  =  0  ,     m  ^  /?. 


Passing  the  constants  outside  of  the  integral  and  dividing  through  by  them  leaves  us  with 

/;^^-f)^^-f)^(i-t)"^=o,  ,„„ 

Taking  the  limit  as  /?  -»  oo,  we  obtain  the  orthogonality  relation  for  the  Laguerre 
polynomials 


(18) 


L(*\x)L(*\x)x:te-xdx  =  Q,     m*n. 


b.     Hermite  Polynomials 

The  Hermite  polynomials,  Hn{x),  are  defined 


(19) 


//„(*)  =  (2*)  Vo 


-«/2,(-if+l)/2         i 


In  a  fashion  similar  to  that  for  the  Laguerre  polynomials,  the  Hermite  polynomials  are 
a  limiting  case  of  the  Jacobi  polynomials  via  the  Gegenbaucr  polynomials.    Specifically, 

nn(x)  =  ni\imrnl2ci"\rll7x) 

/'.-•OO 

which  allows  a  derivation  of  the  orthogonality  relation 


(20) 


oo  2 

—X 


Hm{x)  IIn(x)  e       dx  =  0  ,    m¥=  n 


from  that  of  the  Jacobi  class.  [Ref.  13:  p.  107] 

C.     DISCRETE  EXTENSIONS 

We  turn  now  to  orthogonal  polynomial  classes  which  use  the  discrete  inner  product 
introduced  in  Chapter  II,  Section  A. 2 
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</,*>  =^/W^)w(.v). 

Here  instead  of  a  continuum  the  support  of  the  weight  function  is  concentrated  of  a  fi- 
nite set  of  discrete  mass  points  {  0,  1, ...  ,  A/'}.  The  polynomial  classes  of  particular  in- 
terest are  the  Hahn,  dual  Hahn,  and  the  Racah  polynomials. 
1.     Hahn  Polynomials 

The  Hahn  polynomials  -  actually  discovered  by  Chebyshev  -  were  independently 
realized  by  physicists  working  in  angular  momentum  theory  via  3-j,  or  Clebsch-Gordon 
coefficients  [Ref.  8].   We  define  this  class  by  the  generalized  hypergeometric  series 


(21)  Qn(x',*t  /?,  iV)  =  3/r2 


— n,  —x,  n+a+fi+l 

<x+l,-N 


for  a  >  —  1  ,  /?  >  —  1,  where  N  is  a  positive  integer  and  n  =  0,  \, ...  ,  N.  From  the  power 
series 

we  note  that  the  variable  x  does  not  appear  where  we  have  come  to  expect.   Since 

(-x)k  =  (-x)(-x+l)  ...  (-jc+A-1) 
=  (-l)k(x)(x-l)  ...  (x-k+\), 

we  conclude  that  Qn(x ;  a,  /?,  N)  is  a  polynomial  of  degree  n  in  the  variable  x.  Because 
0  <  n  <  iV,  this  set  of  /V+l  orthogonal  polynomials  is  finite  for  fixed  a  and  /?.  [Ref.  14: 
p.  33] 

The  Hahn  polynomials  satisfy  the  discrete  orthogonality  relation 

v 
(22)  Yj  Qmix  ;  a,  0,  N )  Qn(x  ■  a,  /?,  N  )  w{x  ;  a,  /?,  N )  =  0  ,    m  *  n 

x=0 
where  the  weight  function  is  given  by  a  hypergeometric  distribution 
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fx+a\  [N-x+P 
„(,;*,/*,  A' )  = 

(23) 


(1),  (0,v_x      • 

We  note  that  the  weight  function  may  also  be  written 

{fi+l)N     (<x+l)x(-N)} 


(24)  w(x;a,  /?,  N) 


(l)lV        {l)x(-N-P)x 
We  introduce  the  first  forward  difference  operator  acting  on  x 

Afn(x)=fn(x+l)-fn(x) 
as  a  discrete  analogue  to  differentiation.    Since 


X 


^Afn(t)=fn(x+l)-f„(0), 


[=0 

(a  discrete  analogue  of  the  Fundamental  Theorem  of  Integral  Calculus)  the  first  forward 
difference  operator  is  a  discrete  inverse  of  the  summation  operator.  This  difference  op- 
erator is  used  in  the  Rodrigues'  formula  for  the  Hahn  polynomials,  which  can  also  be 
written  in  terms  of  the  weight  function  w(x;<x,  /?,  N)  with  shifted  parameters  as  was 
done  for  the  Jacobi  polynomials. 

Two  limiting  cases  of  the  Hahn  polynomials  are  of  particular  interest.  In  the 
first  case,  replace  x  by  Nx  in  the  interval  of  orthogonality.  This  in  effect  places  the 
support  of  the  weight  function  on  the  (  N+\)  equally-spaced  points  {  0,  1/A\  2/ A',  ...  ,  1  } 
in  [0,1].  As  N  -*  oo,  the  set  of  discrete  mass  points  tends  towards  the  full  interval 
[0,1]  ;  we  expect  that  this  structure  is  reflected  in  the  discrete  orthogonality  relation  (22) 
becoming  an  integral  orthogonality  with  respect  to  a  continuous  weight  function  on  that 
interval.    Rewriting  the  weight  function  in  (23)  as 

n    ^       r(:c+l+a)      r(N-x+l+P) 

w{x  ;  a,  /?,  A  )  ■■ 


r{x+\)       r(A'-.v+i 

and  using  a  consequence  of  Stirling's  Formula  [Ref.  15:    p.  257] 

55 


Hf+a) 
linn  =  1, 

/-oo  I    (/) 

wc  see  that 

r(a+l)r(/M-l)lim  AT*"* w{Nx ; a,  /?,  /Y)  =  xa  (l-.v/, 

which  we  recognize  as  the  continuous  beta  distribution  for  the  Jacobi  polynomials  nor- 
malized on  [0,1]. 

This  is  indeed  the  case,  and  can  be  easily  verified  by  a  direct  computation  on  the 
hypergeometric  series  definitions  (21)  of  Q„{x  ;  a,  /?,  N)  and  (1)  of  P(*'p)(x)  given  in  Sec- 
tion A.l,  i.e., 

Bm  &(#*;«,  /?,  N)=  — ^—  ft  fi\\-2x). 

Thus  the  Hahn  polynomials  may  be  viewed  as  a  discrete  analogue  and  generalization  of 
the  Jacobi  polynomials.   [Ref.  14:  p.  36] 

The  second  limiting  case  gives  rise  to  an  interesting  class  of  polynomials  which 
has  applications  in  coding  theory  [Ref.  16]. 
a.     Krawtchouk  Polynomials 

For  0  <p<  1,  let  a  =  pt ,  /?  =  (l—p)  t  in  (21),  then  take  the  limit  as  t   tends 
to  infinity  to  obtain  the  Krawtchouk  polynomials,   i.e., 

Ym\Qn(x;pt,{\-P)t,N)  = 


t-*oo 


(25) 


[— n  ,  —x      I    ~| 
-N     ''P   \  =  K»ix'>p>  N) 


[Ref.  14:  p.  38].    A  similar  limiting  argument  applied  to  the  weight  function  shows  that 
the  Krawtchouk.  polynomials  are  orthogonal  with  respect  to  a  binomial  distribution: 

,v 
(26)  Yj  K^x  >P>  N )  K»(x  IP,  N)[     }PX{1  ~pf~X  =  0  ,    m  +  n. 

This  class  of  polynomials  possesses  an  inherent  symmetry,  the  structure  of 
which  can  be  generalized  to  form  other  orthogonal  polynomial  classes.    To  this  end,  we 
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turn  now  to  the  characteristic  of  duality.  For  suitably  defined  functions  u(x)  and  v(.v),  two 
classes  of  orthogonal  polynomials  { pJjMx))  }-.  an^  {  <7„(VM)  j  arc  said  to  be  dual  if 

px(fi(n))  =  qn(v(x)). 

That  is,  interchanging  the  roles  of  the  degree  n  and  the  discrete  variable  x  in  one  class 
produces  the  other.  The  Krawtchouk.  polynomials  provide  an  example  of  an  orthogonal 
class  that  is  self-dual,   i.e., 

Kn(x;p,  N)  =  Kx(n)p,  N). 

This  is  clear  from  the  hypergeometric  series  definition  (25)  given  above.   It  therefore  may 
seem  reasonable  to  suspect  that  there  exists  a  class  of  orthogonal  polynomials  dual  to 
the  Hahn  polynomials.    Such  a  class  does  in  fact  exist,  and  it  is  this  dual  class  which  we 
next  examine.   [Rcf.  17:  p.  657] 
2.     Dual  Hahn  Polynomials 

The  dual  Hahn  polynomials,  Rn(  A(x) ',  y,  d,  N),  are  defined 


(27)  Rn(X(x);V,d,  N)  =  2F2 


—n,  —  x,  x+y+S+l 
V+l,-N 


where     a(x)  =  x(x+y+S+\),     and     n  =  0,  1, ...  ,  N.     These     polynomials     satisfy     the 
orthogonality  relation 


N 


(28)  Yj  Rm(  *M  J  M.  # )  */.(  M*)  \Y.*,N)  M.x  ;  V,S,  N)  =  0,    m  +  n 


x=0 


where 


(29) 


(y+6+l)x((y+S  +  3)l2)x(y+\)x(-N)x(-l)x 
W{*'y'    '      ]  (\)x{(y+d+l)l2)x(d+l)x(y+6+N+2)x 

=        x     (y+S+l)x(y+\)x(-i\')x       2x+y+8+\ 


{\)x(5+\)x{N+y+S+2)x        y+8+\ 


Note  that  Rn(  X(x) ;  y,  5,  N)  is  a  polynomial  of  degree  n   in  the  "variable"  A(x).  The  rea- 
son can  be  seen  directly  from  the  hypergeometric  series 
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K,—  I 

}  [  i-x)  (x+y+S+l)  +j(-x)  +j  (x+y+S+l)  +  /  ] 


By  writing  the  terms  (— jc)*  (.r+y+<5+l)A  in  product  form,  we  obtain 

k-\  k-\ 

\(-x+j)l\(x+y+5+j+l) 

y=o  y=o 

which  becomes 

fc-i 
\(-x+J)(.x+y+d+l+J). 

7=0 

Multiplying  these  factors  as  binomials,  we  have 

k-\ 

n 

y=0 

which  simplifies  to 

k-\ 

n 

y=0 

Taking  only  that  part  which  depends  on  x  yields 

X{x)  =  x{x+y+d+\) 

as  given  above.  [Ref.  18:  p.  48] 

As  discussed  in  the  previous  section,  the  discrete  classes  interact  more  naturally 
with  difference  operators  than  with  differentiation.  To  accomodate  the  quadratic  form 
of  Mx),  we  introduce  the  divided  difference  operator 

,,..,,    a/„(  ;.u»     /„(/  u+ 1  ))-/„(  ;(.vi) 
"^  ','-v,)  =     ^  ;.u,     "       ;.l.r+l,-/(.r,       ' 

The  dual  interplay  between  the  dual  Ilahn  and  Ilahn  polynomials  follows  easily 
from  their  hypergeomctric  definitions: 
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K—  I 

J  [  -jt(x+y+*+l)  +j(y+S+l)  +  j2  ]. 


Rx(A(n);oc,  ft  iV)  =  3F2 


— .v,  — «,  /7+a+/?4-l 

a+l.-.V 


;i 


=  &,(*;  a,  /?,  N). 


This  duality  is  also  reflected  in  the  recurrence  relations  and  diflerence  equations  of  the 
two  classes.    By  interchanging  the  roles  of  n   with  x,    a  with  y,  and  /?  with  d,  the  recur- 
rence relation  for  one  class  leads  to  or  can  be  extracted  from  the  diflerence  equation  for 
the  other  class.   [Rcf.  14:  p.  37] 
3.     Racah  Polynomials 

Both  the  Ilahn  and  dual  Hahn  polynomials  can  be  unified  as  special  cases  of  a 
single  larger  class,  the  standard  notation  of  which  is  similar  to  that  of  the  dual  Hahn 
class.  The  Racah  polynomials,  Rn(  }.{x) ;  a,  P,  y,  S),  are  defined  by  the  hypergeometric 
series 


(30) 


Rn(  ).(x) ;  a,  p,  y,  3)  =  4F3 


'-n,  n+a+p+1,  -x,  x+y+5+l 
a+1,  p+6+l,y+l 


;i 


where  X(x)  =  x(x+y+5+l),  n  =  0,  1, ... ,  N,  and  one  of  a+1,  p+S+l,  or  y+l  equals  —N. 
(Physicists  understand  these  objects  via  6-j  symbols.) 

The  orthogonality  relation  for  the  Racah  polynomials  is 

(31)         Yj  Rm(  A(x) ;  a,  /?,  y,  6)  Rn{  X(x) ;  a,  /?,  y,  S)  w(x  ;  a,  /?,  y,  8)  =  0  ,    m±  n 


where 


x=0 


(32) 


w(x;a,  p,y,8)  = 


(y+6+\)x((y+M)l2)x(a+l)x(p+S+\)x(y+l)x 

(i)x{(y+d+i)l2)x(y+d-a+i)x{y-p+i)x(6+\)x 

(y+5+l )x  (a+ 1  )x  (P+S+ 1  )x  (y+ 1  )x     2x+y+6+ 1 


(l)x(y+S-a+l)x(y-P+l)x(d+l 


y+d+l 


[Ref.  19:  p.  24].  Note  that  w(x  ;  a,  P,  y,  <5)  has  a  "well-poised"  structure.  This  means  that 
the  pairwise  sum  of  numerator  and  denominator  parameters  is  constant,  i.e., 
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(y+c5+l)  +  (l)  =  y+c5+2 


y+<5+3 


+ 


y+<5+] 


=  y+S+2 


(a+1)  +  (y+d-a+l)  =  y+S+2 

((J+S+\)  +  (y-P+l)  =  y+S+2 

{y+l)  +  {S+l)  =  y+S+2. 

The  same  could  be  said  for  the  weight  function  of  the  dual  Hahn  polynomials.  (To  say 
that  a  hypergeometric  series  is  well-poised  means  that  for 


r+\^r 


a,, ...  ,  ar+1 
_  bx, ...  ,br 


\x 


the  parameters  a  and  b   satisfy  the  following  relation 

a,4-l  =  a2+bx  =  •••  =  ar+l+br 

Well-poisedness  is  an  important  property  of  certain  summable  hypergeometric  series.) 

From  the  hypergeometric  series  definition  (30)  of  the  Racah  polynomials,  we 
find  that  when  the  roles  of  x  and  n,  a  and  y,  and  /?  and  S  are  all  interchanged,  the  series 
is  unchanged.  Thus,  like  the  Krawtchouk  polynomials,  the  Racah  polynomials  are  self- 
dual. 

To  recover  the  Hahn  polynomials  as  a  limiting  case  of  the  Racah  polynomials, 
let  y+l  =  —N  and  5  -*  oo.  Thus  by  formulas  (30)  and  (21), 


Jim  4F3 


-n,  n+a+p+\,  -x,  x+y+S+\ 
a+1,  fi+S+l,  -N 


;i 


3^2 


— n,  n+a+fi+l,  —x 

a+1,  -N 


;l 


and  then  by  formulas  (32)  and  (24), 


(Da 


\\mw{x;a,  (J,y,S)  = 

<5-»oo  IP+U,V 


w 


(*;a,  ft,  N). 


Likewise,  by  letting  a+1  =  —N  and  P~*  oo,  we  obtain  the  dual  Hahn  polynomials. 

Figures  5  and  6  together  with  Tables  2  and  3  provide  a  hierarchy  of  classes  dis- 
cussed in  this  chapter.  Tables  4-15  summarize  information  about  selected  classes  in  this 
hierarchy  [Refs.  14,  12,  20,  21  ,5,  13,  9  ,  22]. 
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Racab 


Laguerr 


Hermtte 


Dual  Hahn 


Meixner 


Figure  5.      Hierarchy  of  Limiting  Cases 
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Jacobi 


Chebyshev     Chebyshev     Legendre 
First  Kind    Second  Kind 


Figure  6. 


Hierarchy  of  Special  Cases 
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Table  2.     TABLE  OF  CLASSES 

Class 

Symbol 

Ilypergcometric  Scries 

Racah 

Rn(  A(x)  ;  a,  fi,  y,  <5) 

r  V-n,  n  +  a  +  p  +  1 ,    -x.x  +  y  +  5  +  I      ' 
4^|_              a+l,P  +  S+[,y+l                'j 

Dual  Hahn 

Rn(A(x);y,S,  N) 

r  [~— n,  —x,  x  +  a  +  P  +  1     . 
2t2\_          y+  1,    -N           >l_ 

Hahn 

Qn(x;a,  p,  N) 

r  [~— x,    —n,  n  +  a  +  P  +  1     . 
^\           a  +  1,    -iV            '*_ 

Jacobi 

pf%) 

(a+l)„      rf-A7,  /?  +  a  +  /?  +  1       l-.r" 
«!         2   ][_            a+1             '2 

Krawtchouk 

Kn{x\p,  N) 

2*i|_    "l^v  X'~\ 

Meixner 

Mn(x ;  p,  c) 

2f'L  %  Y; l~~J 

Laguerre 

t(x)(  ^ 

n!         '    ' 

_a+"l;j:] 

Charlier 

Cn(x  ;  a) 

2y  o|_       _        »       a  J 

Hermite 

Hn(x) 

^r2Fo[-nl2A-n+W._y 
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Table  3.     TABLE  OF  CLASSES  (CONTINUED) 


Gegenbauer 

Cn  (x) 

W)„     r[-n,  n  +  2X      \-x 

n\      2   1    X  +  ll2     '       2      J 

Chebyshev, 
First  Kind 

Ux) 

rY-n,  n.    \  -x  ~] 
2tll   1/2    '       2      J 

Chebyshev, 
Second  Kind 

Un(x) 

("+l)a*iL      3/2       ;       2      J 

Legendre 

PnW 

r  f-rt,    «  +  1  .      1  -  X 

2t]l        1         '       2      J 
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Table  4.     RACAH  POLYNOMIALS 


Symbol:    Rn(  X(x)  ;  a,  /i,  y,  S) 

Interval:    For  N  a  positive  integer,  x  =  0,  1, ... ,  N. 

Weight: 

w(x  ;  a,  /?,  y,  (3)  = 


(y+<5+lU(y+<5  +  3)/2)x(a+l)x(/?+<HlUy+l)x 


(1),  ((y+o+\)l2)x  (y+d-a+l)x  (y-P+\)x  (d+\)x 


Norm: 


£  [  /*„(  iW  ;  a,  fi,  y,  5)  f  w(x  ;  a,  /?,  y,  <5) 


x=0 

«!  („+«+/?+!)„  ( /?+!)„  (a-<5+l)„  (a+/?-y+l)„ 


where  if 


or  if 

0+5+1  =-N, 


or  if 

y+l  =  -iV, 


{a+P+2)2n{(x+\)n(P+d+i)n{v+\)n 


(y+6+2)N(-P)N 
1    '"'  (7-P+Dn(M)h  ' 


(y+<5+2).v((5-a),v 

M  = ■ ■ — 

{y+d-<x+l)ff{6+[)N  ' 


M);V(«+/?+2),v 
=  (a-<5+l)„  (/?+!)„ 
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Table  5.     RACAH  POLYNOMIALS  (CONTINUED) 


Hypergcometric  Scries: 

niu  \  a        s\        rY—n,n+a+P+\,—x,x+y+8  +  \     ,~1 

where  X{x)  =  x(x+y+S+\),  n  =  0,  1, ...  ,  N, 
and  one  of  oc+1,  fi+d+1,  or  y+1  equals  —  N. 

Recurrence  Relation: 

A(x)  Rn(  X{x) ;  a,  /?,  y,  8)  =  a(«)  [  /?n+1(  X{x) ;  a,  /?,  y,  d)-Rn{  X{x) ;  a,  0,  y,  8)  ] 

-c(«)  [  RJL  X{x) ;  a,  p,  y,  S)-Rn_,(  X(x) ;  a,  fi,  y,  8)  ] 

where 

(n+a+P+\)  (n+a+1)  jn+p+M)  (n+y+l) 
a{"}  ~  (2n+a+P+l)  (2n+a+P+2) 

n  (n+p)  jn+a+P-y)  (n+cc-S) 
C(/7)_        (2n+<x+P+l){2n+a+P) 

Difference  Equation: 

;.(«)  Rx(  X(n)  ■  a,  p,  y,  8)  =  A(x)  [  Rx+X(  X{n) ;  a,  p,  y,  8)-Rx(  X(n) ;  a,  p,  y,  8)  ] 

-C(x)  [  Rx(  A{n)  ;  a,  /?,  y,  6)-Rx_x(  X{n)  ;  a,  /?,  y,  5)  ] 


where 


(.v+y+d+l)  (-v+y+1)  (x+8+p+l)  (x+a+1) 

W  "  (2jc+y+<5+l)  (2*+y+<5+2) 

_  g  (.t+<5)  (x+y+O-oc)  (x+y-/?) 

W  (2x+y+<5+l)(2x+y+()) 
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Table  6.     DUAL  HAHN  POLYNOMIALS 


Symbol:    K„(  X(x) ;  y,  5,  N) 

Interval:    For  N  a  positive  integer,  x  =  0,  1, ...  ,  N. 

Weight: 

c     _       (y+S+\)x((y+6+3)j2)x(y+\)x(-N)A-l)X 
wlx  '  V   o     /V  )  = 

1    'A    '      )          (l)^((y+5+I)/2)x(5+l)JC(y+<5+iV+2)JC 

Norm: 

/V                                                                                                                                 _ 

£[*„U(;c);y,<5,  .V)]2Htv;y,c5,  tf) - 

x=0 

\  n  ) \   N-n   ) 

-l 

Hypergeometric  Series: 

Rnmx);r,S,N)  =  3F2[-"'-£*+j£S+i--1] 

where  yi(jr)  =  x  (x+y+S+l) 

Recurrence  Relation: 

-A{x)  Rn(  A(x)  ;  y,  <5,  AT)  =  Z?(*)  [  K„+1(  ;.(-r)  ;  y,  <5,  N)-Rn(  X(x)  ;  y,  <5,  N)  ] 

-Z)(/?)  [  /?„(  ;.(*)  ;  y,  (5,  iV )-*„_,(  /(x) ;  y,  0,  tf)  ] 

where 

Z?(«)  =  (JV-«)(y+l+/*) 

£>(/?)  =«(  N+l+S-n) 

Difference  Equation: 

-nRn{  ).{x)  ;  y,  <5,  iV)  =  &(*)  [  K„(  J(x+1)  ;  y,  5,  iY)-^(  A(*)  ;  y,  5,  N)] 

-d{x)  [  Rn{  k{x)  ;  y,  <5,  N  )-Rn(  A(x-l)  ;  y,  S,  N  )  ] 

where 

{x+y+S+l)  (x+y+ 1 )  (  N-x) 

(2x+y+S+l){2x+y+S+2) 

x  (x+5)  (x+y+5+N+ 1 ) 
(2x+y+S)  (2x+y+d+\) 
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Table  7.     HAHN  POLYNOMIALS 


Symbol:   Q„{x  ;  a,  /?,  N) 

Interval:    For  N  a  positive  integer,  .v  =  0,  1, ...  ,  N. 

Weight: 

t  a    *n      (x+<x\(N-x+0\      (/?+D,y     (a+l)x(-AQ; 


Norm: 
00 
£[&,(*;  a,  /?,  iY)]2wU;a,  /?,  N) 

x=0 

a+/?+ 1  \        N        )\  n 


2/7+a+/?+l  AV\ 

r(a+i)r(a+/?+i)       r(rt+/?+i)n«+i) 


r(/?-hl)  r(«+a+l)r(«+a+yS+l) 


Hypergeometric  Series:    For  a,  /?>  — 1  and 
«  =  0,  I,...  ,/V, 


fl^;«,AiV)-3^["n'^^'"*J1} 


X 
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Table  8.     HAHN  POLYNOMIALS  (CONTINUED) 


Recurrence  Relation: 

-xQn(x;a,  P,  N)  =  b(n)  [  Qn+l(x  ;  a,  p,  N)-Qn(x;a,  P,  N)} 

-d{n)  [  Qn(x  ;  a,  p,  N )  -  g^fr  ;*,  P,  N)] 

where 

(w+a+/?+ 1)  (n+a+1)  (  N-n) 
W"     (2«+a+/?+l)  (2«+a+0+2) 

n(n+p)  (n+a+P+N+\) 
{n}"  (2n+o.+p)(2n+a+p+l) 

DifTerence  Equation: 

-).{n)  Qn(x  -a,  P,  N)  =  B{x)  [  Qn{x+ 1  ;  a,  0,  N )  -  &(*  I  «,  /?,  # )  ] 

-Z)(.v)  [  &(* ;  «,  /?,  iV)-a(x-l;a,  /?,  iV  )  ] 

where 

5W  =  (Ar-Jc)(a+l+Jc) 

Z)(jc)  =  jc(/Y+1+/?-jc) 

A(/?)  =  «(«+a+/?+l) 

Rodrigues'  Formula: 

wrao^*" 

=(^)A„ 
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Table  9.     JACOBI  POLYNOMIALS 


Symbol:  /*•  '\x) 
Interval:  [-1,1] 
Weight:   (1-*)"(1  +  xf 

n  +  a\       (<*  +  1)„ 


Standard 
Norm: 

•J-i 


ization:   JPl'/?,(l)  =  ^    „    J= ^p 


r/^  %)?  (i  -  *r  (I  +  xf  dx  =  ^+1r(,+a+i)r(,+/?+i) 

1/1       WJ   U        J   l    +    ;  (2A/+a+/?+l)A?!r(«+a+/J+lj 


Hypergeometnc  Series:    /*'(*)  = — - — 2FA  a+1  ;  ~T~ 

Recurrence  Relation:   2  («+l)  («+a+/?+l)  (2«+a+jff)  /*;/V) 

=  (2«+a+/?+l)  [  (2«+a+/?)  (2/i+a+/?+2)  jc+a2-/?2  ]  Pj^Gr) 

-2  («+a)  (*+/?)  (2/7+a+/?+2)  Pj^Cx) 

Differential  Equation: 

(l-*2)  /  +[/?-«-  (oc+0+2)  *]/+«  («+a+/?+l).y  =  0  ,  y  =  P<*,/?)(*) 

Rodrigues'  Formula: 

2"  n\  P%>  p\x)  =  (- 1 )"  ( 1  -x)-*  ( 1  +*)"*  -~p  [( 1  -  *)a+"(  1  +  xf+n] 

dx 

Generating  Function: 


V  /f  ^V)  r",  2  ([-t+Jl-2xt  +  t2)-*  (i+r+v/l-2-v/  +  /2V/? 

n-o  v/(1-2av+/2)     V  7      V 

Explicit  Expression:   /*"'  "(*)  =  2""  V  ("^a)  (^)  (jc-1)-*  (x+1)* 
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Table  10.     LAGUERRE  POLYNOMIALS 


Symbol:    L£\x) 
Interval:   [0,  oo) 
Weight:   x*  e~x ,     a  >  —  1 


i-\Y 


-rr^^^.-^.I&^+O 


Standardization:    L„  (1)  = : — x"+ 

Norm: 

Hypergeometric  Series: 

Recurrence  Relation: 

(if+1)  LJJ,(*)  =  [  (  2n+a+l  )-*  ]  L^W  -  («+«)  *&(*) 

DifTerential  Equation:   jcj/"  +(  a+l—  .v)  j/+/?j/  =  0  ,  y  —  Ll*\x) 
Rodrigues'  Formula: 

n\  x   e         ax 
Generating  Function: 

CO 

^L<»>W/"  =  (l-rr-'  exp(-^j-) 


n=0 


Explicit  Expression: 


tfw-I>»*  (£*)£** 


*=o 


71 


Table   11.     HERM1TE  POLYNOMIALS 

Symbol:    H„{x) 

Interval:   (—00,  00) 

Weight:   e~*2 

Standardization:    //„(!)  = 

=  2  x  +  ... 

Norm: 

[  Hn(x)  f  e~x  dx  =  Jn  2n  n\ 

Hypergeometric  Series: 

IIn(x)  =  (2xf2F0 

"— «/2  ,  (-/i+l)/2.        1 

-              '        2 

Recurrence  Relation: 

Hn+1(x)  =  2xHn(x)-2nHn_](x) 

DifTerential  Equation: 

y"-2xy'  +  2ny  =  0,  y=Hn{x) 

Rodrigues'  Formula: 

dx     u          J 

Generating  Function: 

V    HiM'         ix,-,2 
L        „!        — 

«=0 

Explicit  Expression: 

fe=0                      v            y 

72 


Table  12.     GEGENBAUER  (ULTRASPHERICAL)  POLYNOMIALS 


Symbol:    Q<%v)  (or  P?\x)) ,    ^  > -1/2 

Interval:    [-1,1] 

Weight:   (1  -  x2)^'2 

Standardization: 

JS  cf(i)  -^W^-0'1-2- 

Norm: 

% 

['[crw]^,-^-'^-2?'^;:;-"''2:12 

_!  L    '       J    v           ;                       (n+X)  n\  r(«+2>.) 

Hypergeometric  Series: 

r%  .^      (2y^«     rf-«,/i+2A.    l-x 
n  {  >~     n\      2   '  [     i+1/2     '      2 

Recurrence  Relation: 

(//  +  1)  ($,(*)  =  2  («+/l)  *  C<;)(.v)  -  (n+ll-1)  c^,w 

DifTerential  Equation: 

(l-x2)y"-{2A+l)xy'+n(n+2A)y  =  0,  y=C(n%x) 

Rodrigues'  Formula:    C0(;,(.v)  =  1  ,    C/V)  =  2Xx, 

r  n\  (i+1/2),  (l-.t2);-,/2  C^Xx)  =  (-If  (2X)n  •£-  [  (1-x2)^-'/2] 

Generating  Function: 

oo 

^crw,"=(i-2«+,2)-;- 

n=0 
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Table  13.     CHEBYSHEV  POLYNOMIALS  OF  THE  FIRST  KIND 


Symbol:    T„(x) 

Interval:   [-1,1] 

Weight:   (1  -.v2)-1/2 

Standardization:    r„(l)=l 

Norm: 

*1 

r 

i. 

I 

[7-„W]2(i-,r"2*={f^fo 


Hypergeometric  Series: 

UD-Jf,[7/i":iTi] 

Recurrence  Relation:    Tn+l(x)  =  2xTn(x)  —  r„_,(jr) 
Differential  Equation:   (1  —  x2)y"  —  xy'  +  n2y  =  0  ,     y  =  Tn(x) 
Rodrigues'  Formula: 

Generating  Function: 

oo 

yr„Wf"=     1-*f  2  ,  -i<*<i,  |/|<i 

^  l-2.v/+f2 


n=0 


Explicit  Expression: 

Tn(x)  =  cos(nO)    with   jc  =  cos  0  ,  0  ^  0  ^  tt 


t^       n   V  ^     ,^  («-/r-  1)!  „_2/< 


k=0 
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Table  14.     CHEBYSHEV  POLYNOMIALS  OF  THE  SECOND  KIND 

Symbol:    U„(x) 

Interval:    [—1,1] 
Weight:   (1  -  x2)1'2 
Standardization:    Un(l)  =  n+l 
Norm: 

-i 

Hypergeomctric  Series: 

C/„W  =  (/t+l)2F1[-rt^+1;-^] 

Recurrence  Relation:    Un^(x)  =  2xUn(x)  —  Un_{{x) 

Differential  Equation:   (1  —  x2)y"  —  3xy'  +  n  (n+2)y  =  0  ,  y  =  Un(x) 

Rodrigues'  Formula: 

Generating  Function: 

Yun(x)tn  = l- =-,    -1<jc<1,     U|<1 

to  1-2  W 


Explicit  Expression: 


sin(«  +  1)0 

UJx)  = : — with    x  =  cosd  ,  0<d  <n 

sin  0 

i<=o 
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Table  15.     LEGENDRE  (SPHERICAL)  POLYNOMIALS 


Symbol:    Pn{x) 

Interval:    [-1,1] 

Weight:    1 

Standardization:    P„(l)  =  1 

9 

.\orm.   J1^  lPn(x)j  «-*  -  2n+\ 

Hypergeometric  Series:    Pn{x)  = 

,,[- 

-«,  «+l      1— X 
1        '      2     J 

Recurrence  Relation:   (n+1)  Pn+ 

w- 

(2/7+1)  xPJtf-nP* 

-iM 

Differential  Equation:   (1— jc2)j, 

'"  -  2. 

xy'  +  n(n+l)y  =  0, 

y  = 

Pn(x) 

Rodrigues'  Formula:   P„(x)  =  — 

-ir 

i*  [<i-*T] 

Generating  Function: 

«=0 

-2xrt 

-/2)-"2,    -1<*<1 

,   1 

t\<\ 

Explicit  Expression:    Pn(x)  =  -rr 

[«/2] 

Z( 

*=0 

-^©(W 

2* 
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V.     APPLICATIONS 

Orthogonal  polynomials  and  special  functions  in  general  have  been  studied  exten- 
sively in  mathematics  and  other  fields  since  the  eighteenth  century.  Presented  below  are 
a  few  of  the  traditional  applications  of  selected  classes  of  orthogonal  polynomials.  Our 
first  few  applications  come  from  numerical  analysis. 

A.     ECONOMIZATION  OF  POWER  SERIES 

Economization  of  power  series  is  a  technique  used  to  reduce  the  degree  of  a 
polynomial  approximation  to  a  given  function. 

The  maximum  norm  (or  L°°  —  norm)  for  a  continuous  function  on  a  compact  inter- 
val [a,b~\  is  defined  as 

ll/IL  =  max  |/(.v)|. 

xe  [a,b] 

This  norm  is  not  induced  by  an  inner  product,  but  nevertheless  has  many  uses  in  applied 
mathematics,  including  numerical  analysis. 

The  minimax  property  of  the  Chebyshev  polynomials  states  that  of  all  n"'  degree 
monic  polynomials  (i.e.,  leading  coefficient  1),  2*-"Tn(x)  has  the  smallest  maximum  norm 
on  [  —1,1]  [Ref.  23:  p.  106].  The  justification  for  this  statement  is  deferred  until  Section 
C.  Hence  the  best  approximation  in  the  maximum  norm  to  the  function  x"  on  [  — 1,1] 
by  a  function  of  lower  degree  is  f„(x)  =  xn  —  2l~"TJx),  So,  given  a  function  and  a 
polynomial  approximation  to  that  function  (e.g.,  from  a  Taylor  series  expansion),  suc- 
cessively replace  the  highest  powers  x"  mthfn(x)  to  obtain  a  polynomial  approximation 
of  lower  degree  [Ref.  23:  p.  125]. 

Example  1:    Let/(.xr)  =  sin  x.  The  Maclaurin  series  for  this  function  is 


X  =  L  T? 


(-ir     2*4-1 

sin  x  =   7     — — ——  x 


(2k  +  1)! 

and  this  series  is  convergent  for  x  e  U.  If  truncated  after  the  x5  term,  the  polynomial 
approximation  for  the  function  is 

3  5 

smxxx-  4r  +  ^r 
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with  a  maximum  error  of  0.0002  for  jc  e  [  — 1,1].  We  use  the  fifth  degree  Chebyshev 
polynomial  T5{x)  =  L&c5  —  20.y3  +  5x  to  obtain 


5  r     3  c 

X  DX  JX 


5!  480  1920 

with  an  error  not  exceeding 

max  I  T5(x)  I  i 

in  [  — 1,1].  Thus  the  approximation 

_    *L   i_fj*l  5.v    \       1915      _  _75_    3 

smx~x       3!    +^480  1920  )       1920  ^       480  X 

has  an  error  whose  magnitude  in  [  — 1,1]  does  not  exceed 

0.00052  +  0.0002  =  0.00072. 

Compare  this  with  the  maximum  error  of  0.00833  for  the  Maclaurin  series  which  is 
truncated  after  the  x*  term.  For  a  cubic  polynomial  approximation  of  sin  x,  the  "econ- 
omized" polynomial  has  a  maximum  error  that  is  significantly  smaller  (less  than  one 
tenth)  than  that  of  the  truncated  Maclaurin  series. 

The  next  three  applications  illustrate  the  usefulness  of  the  zeros  of  orthogonal 
polynomials.  The  first  two  come  from  numerical  analysis,  the  third  from  a  problem  in 
electrostatics.  We  begin  with  a  preliminary  discussion  of  a  fundamental  technique  from 
numerical  analysis. 

B.     POLYNOMIAL  INTERPOLATION 

Polynomial  interpolation  is  a  method  of  approximating  a  given  function  with  a 
polynomial  that  matches  (interpolates  )  the  function  at  specified  points  (called  nodes  or 
abscissae)  x{,  ...  ,  xn  [Ref.  24:  p.  497].  Given  a  function /(.v)  and  n  distinct  nodes  in  a 
compact  interval  [«,/>],  there  is  a  unique  polynomial  of  degree  (n  —  1)  that  passes 
through  the  points  (.xr, ,  fix,)),  1  <  /  <  n. 

For  each  /  =  1,2, ...  ,  n,  define  a  polynomial  of  degree  (n  —  1)  by 
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n 


x- 


Xi—Xi 

k±i 

Clearly,  for  each  k  =  1,  2, ...  ,  n,  ntJxk)  =  0  \f  k¥*  i  and  rc^fo)  =  1,  that  is,  nttn(xk)  =  bik. 

n 

(We  may  thus  equivalently  express  nlin(x)  =  ln(x)l({x—x,)  l„'{x,)),  where  ln(x)  =  n  (-v— -ya)  is 
the  unique  monic  polynomial  of  degree  n  that  vanishes  simply  at  each  node  xh, 
*-l,  2, ...,«.) 

The  Lagrange  interpolating  polynomial   is  now  given  by 


Clearly,  L{(jc*)  =/(.yJ  for  &=  1,  2, ...  ,  n,  and  uniqueness  of  the  degree  (n  —  1)  interpolat- 
ing polynomial  is  guaranteed  by  the  Fundamental  Theorem  of  Algebra. 

C.     OPTIMAL  NODES 

In  this  section  we  address  the  issue  of  estimating  the  maximum  size  of  the  interpo- 
lating error  \\f—  L{|L.  Assuming  /  is  suitably  differentiable  in  \_a,b~\,  it  can  be  shown 
[Ref.  25:  p.  188]  that  there  exists  a  value  Kn  (which  depends  on/^'in  [a,6])  such  that  for 
any  x  in  [a, 6], 

\f{x)-Lfn{x)\<-^-  \ln{x)\, 

where  ln(x)  =  Yl(x  —  x,).  (Note  again  that  ln(x)  is  a  monic  polynomial  of  degree  n  which 
vanishes  at  the  nodes.)  This  implies  that  \\f—  l/Jl^  is  minimized  by  making  the  optimal 
choice  of  nodes  xlt ...  ,  x„  in  [a, b~\  which  minimizes  ||  /„  H^.  Surprisingly  perhaps,  this  hap- 
pens precisely  at  the  zeros  of  the  Chebyshev  polynomials  Tn{x),  scaled  to  the  interval 
[a,/?].  We  sketch  the  reasons  below. 

For  simplicity,  we  take  our  interval  of  interest  to  be  [  — 1,1]  instead  of  ia,b~]  without 
loss  of  generality.    The  transformation 

*=,i=2.\  t+fb+a 


or  equivalently, 
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'-2te)-'- 


is  a  one-to-one  continuous  mapping  between  the  intervals  [  —1,1]  and  \_a,b~\. 

Now  let  xi  be  such  that  T„(x,)  =  0,  /  =  1,  2, ...  ,  n  (see  Chapter  I,  Section  A. 3),  i.e.,  let 
L(x)  =  ~^T  Tn{x).  From  the  recurrence  relation  (Equation  (4)  in  Chapter  I),  Tn(x)  has  a 
leading  coefficient  of  2"_1  ;  hence  this  ln(x)  is  monic.  Moreover,  from  the  definition 
Tn{  cos  0)  =  cos  nd,  it  follows  that 

(i)  tw-c-rf-jb- 

where  yt  =  ytJ,  =  — z —  ,  /  =  0,  1, ...  ,  n.  Now  suppose  mn{x)  is  another  monic  polynomial 
of  degree  n   such  that 


(2)  II  mn  i^  <  ||  /„  IL  =  ^rr  • 

Combining  (1)  and  (2)  we  see  that  we  must  have  for  /  =  0,  1, ...  ,n, 

mn(yi)  <  tniy'i)  =  — ~  if ' is  even> and 

0)  2_! 

mn(yi)  >  lM  =  -^~r  if  * is  odd- 

Thus,  the  polynomial  pn(x)  =  mn{x)  —  l„{x)  has  degree  at  most  (n  —  1)  (since  it  is  a  differ- 
ence of  two  monic  polynomials  of  degree  n)  with  at  least  n  zeros,  one  in  each  interval 
(y,  ,yi+l),  0  <  /<  n  —  1,  by  (3).  This  contradicts  the  Fundamental  Theorem  of  Algebra, 
and  so  no  such  polynomial  mn(x)  satisfying  (2)  exists.  Hence  the  choice  of  nodes  xt  de- 
termined by  ln{x)  =  -r— -  T„(x)  minimizes  ||  /  \\,x  =  max  |  ln(x)  \  over  all  possible  monic 
polynomials  ln(x),  and  therefore  over  all  possible  choices  of  interpolating  nodes  x,. 

It  should  be  emphasized  that  this  Chebyshev  interpolation  allows  an  a  priori  error 
bound  for  all  x,  but  is  not  always  best  possible  for  every  x  using  other  interpolation 
schemes.  (For  example,  if  equally  spaced  nodes  xif ...  ,xn  are  used,  then  trivially,  the  er- 
ror/(jc,)  —  LJn{x)  =  0,  even  if  x  =  x,  is  not  a  Chebyshev  zero.)  Chebyshev  zeros  are  opti- 
mal when  one  has  freedom  in  the  choice  of  nodes.  For  details  regarding  the  practicality 
of  Chebyshev  interpolation  and  some  asymptotic  results,  see  [Kef.  24]. 
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D.     GAUSSIAN  QUADRATURE 

Quadrature  formulas  are  used  in  that  area  of  numerical  analysis  concerned  with  the 
approximate  integration  of  a  function /(.v)  against  a  weight  function  w(x)  >  0  on  an  in- 
terval (a,b),  when  the  explicit  evaluation  is  intractable. 

An  interpolator);  quadrature  is  such  a  rule  that  uses  interpolating  polynomials,  such 
as  Lagrange  polynomials: 

f  /(*)  w(x)  dx  *  f    lfn{x)  w(x)  dx, 

n 

where  Lfn{x)  =  X/(-*"<)  *,(*),  as  described  above.   This  can  be  rewritten  as 

n 

(1)  I   f{x)w(x)dx^yf(xl)wh 
where  the  weights  wt  are  given  by 

(2)  wt  =       nln{x)  w(x)  dx. 

Thus,  interpolatory  quadrature  is  basically  a  weighted  sum  of  the  function  values /(x,) 
at  the  nodes  xt,  i=  1,  2, ...  ,n,  as  are  numerical  integration  recipes  such  as  Simpson's 
Rule  and  the  Trapezoidal  Rule. 

For  specified  nodes  xu  ...  ,  xn  e  \_a,b~\,  the  n  weights  w,, ...  ,  wn  computed  in  (2)  for  the 
quadrature  (1)  will  be  exact  for  polynomials  fe  P„-l[a,b'],  but  we  can  do  better.  [Ref. 
25:  p.  236] 

In  Gaussian  quadrature,   we  ask  for  the  location  of  the  n  nodes  xu  ...  ,  x„  as  well  as 
the  n  weights  w,, ...  ,  wn  in  order  for  the  quadrature  rule  (1)  to  be  exact  for  polynomials 
/e  ^2n-iCa^]-  At  first  glance,  this  seems  to  be  an  extremely  complicated  computational 
problem,  but  the  solution  falls  out  simply  when  the  theory  of  orthogonal  polynomials 
is  applied. 

Let  \pn{x)}Zn  be  the  class  of  polynomials  orthogonal  with  respect  to  the  weight 
function  vv(.v)  on  [_a,b~],  say,  by  Gram-Schmidt,  and  let  Lfn{x)  be  the  Lagrange  polynomial 
that  interpolates  fe  P^-iD*^]  at  l^ie  zeros  of  pn(x),  so  that/(.v)  -  Lf„(x)  e  /\n_,[a,/>]  van- 
ishes at  the  zeros  also.  Since  the  zeros  of  pn{x)  are  real,  simple,  and  lie  in  (a,b),  we  have 
the  property  that 
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f{x)  -  Lfn{x)  =  pn{x)  qn_x{x), 
where  qn_x  e  Pn-X\_a,b~\,  Hence,  by  orthogonality, 

f{x)  w(x)  dx  -       Lfn{x)  w{x)  dx  =       pn(x)  qn_x{x)  w(x)  dx  =  0. 
Ja  Ja  Ja 

We  have  thus  shown  that  the  resulting  n-point  Gaussian  quadrature  rule  is  exact  for 
/s  Pin-^b].  Moreover,  by  a  theorem  of  Stieltjes,  if/(jc)  is  continuous  on  a  finite  inter- 
val [a,&],  then 


I"  /  f 

lim      LUx)  w(x)  dx  =      fix)  w(x)  dx 


By  applying  the  ChristofTel-Darboux  formula  (Equation  (9)  of  Chapter  III  with 
notation  from  Chapter  III,  Section  B)  and  using  (1),  it  is  possible  to  derive  an  alternate 
expression  for  the  weights  (2): 


"nPn'(xi)Pn+M'i) 

which  are  referred  to  as  the  Christoffel  numbers.  [Ref.  8  ] 

Thus  Gauss-Jacobi,  Gauss-Chebyshev,  and  Gauss-Legendre  are  the  names  given  to 
Gaussian  quadratures  involving  the  weight  functions  and  orthogonal  polynomials  from 
the  Jacobi,  Chebyshev,  and  Legcndre  classes,  respectively. 

E.     ELECTROSTATICS 

The  zeros  of  the  Jacobi  polynomials  play  an  interesting  role  in  a  problem  of  Stieltjes 
concerning  electrostatic  equilibrium.  In  this  problem,  fix  "masses"  of  positive  charge  a 
and  P  at  the  points  x  =  1  and  x=  —  1,  respectively.  Then  place  n  point  masses  of  positive 
unit  charge  in  the  interval  (  — 1,1)  so  that  they  are  free  to  move.  These  interior  masses 
are  now  subject  to  a  "logarithmic  potential",  that  is,  a  repelling  force  that  is  proportional 
to  the  logarithm  of  the  distance  separating  them.  The  problem  is  to  determine  the  dis- 
tribution of  the  point  masses  x, ,  /  =  1,  2,  ...  ,  n,  when  the  system  is  in  equilibrium. 
Mathematically,  this  is  equivalent  to  maximizing  the  force  function 


82 


(3)  F(x„  ...  ,  xn)  =  a  £  log(l-^)  +  /?  V  log(l+^)  +^]      log  I  ^-.r;  | . 

i=  1  (=1  1  <i<j<n 

The  logarithmic  terms  give  the  restrictions  xx*fc  —  1,  x„¥=  I,  and  xt=fcXj  for  /#,/.  This 
otherwise  continuous  function  gives  the  equilibrium  points  by  setting 

(4)  4^  =  0,  Z-1,2,...,/! 

ex) 

and  solving  for  xlf ...  ,x„.  To  solve  this  system  of  a;  nonlinear  equations  in  n  unknowns, 
Stieltjes  introduced  the  polynomial 

n 
Pn(x)  =  ||  (x~xi) 

and  reduced  (4)  to 


/=i 


1  Pn   (*l)  P  a  r»       •        1    O 

~ 77—-  +  -7- : =  0,  1=  1,2, ...  ,  « 

2  pn'{Xi)  l+Xi         i-xl 


which  becomes 

2\    . 


(l -*>„-(*,)  +  [  2/?  -  2a  -  (2a  +  2/?)  jr,  ]  />„'(*/)  =  0 

for  /=  1,2, ...  ,  n. 

Since  the  polynomial 

(l-x2)  pn"(x)  +  [  2/?  -  2a  -  (2a  +  2/?)  x  ]  ,„'(*) 

is  of  degree  at  most  n  and  vanishes  at  x  —  xt,  i=  1,  2, ...  ,  n,  it  can  be  set  equal  to  a  scalar 
multiple  Xpn{x)  which  also  vanishes  at  these  points.    Hence 

( 1  -  x2)  pn"(x)  +  [  2/?  -  2a  -  (2a  +  2/?)  *  ]  Pn'{x)  -  ^„(jc)  =  0. 

Attempting  a  power  scries  solution  to  this  second-order  differential  equation  leads  to  the 
observation  that  polynomial  (i.e.,  terminating)  solutions  exist  if  and  only  if 

X  =  -n(n  +  2a  4- 2/?-  1). 
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Rearranging  terms  yields 

(1  -  *2)  Pn'\x)  +  [  (2P  -  1)  -  (2a  -  1)  -  ((2a  -  1)  +  (2/?  -  1 )  +  2)  x  ]  pn'(x) 
+n(n  +  (2a  -  1)  +  (2/?  -  1)  +  1)  pn{x)  =  0 

which  is  the  differential  equation  for  Ff'~l,2p~l\x).  Thus  the  equilibrium  positions  of  the 
unit  charges  occur  at  the  zeros  of  this  Jacobi  polynomial.   [Refs.  8,  13:  p.    140] 

The  zeros  of  the  Laguerre  and  Ilermite  polynomials  can  be  developed  as  solutions 
to  similar  electrostatic  equilibrium  problems. 

F.     SPHERICAL  HARMONICS 

We  investigate  another  application  of  orthogonal  polynomials  to  problems  in 
mathematical  physics.  In  Cartesian  coordinates  (x,y,z),  the  Laplacian  operator  of  a 
function  u(x,y,  z)  is  defined  as 

a2  ,2  -,2 

^2  G  U      ,       O  U      ,      O  U 

V  u  =  — r  +  — r  + 


dx  dy  dz 


When  converted  to  spherical  coordinates  (r,  0,  </>),  this  operator  acting  on  a  function 
u(r,  6,  0)  becomes 


(5) 


V  u  =  — : 


(r2  w,  )r  + 


1 


sin  0 


( sin  6  u6)e  + 


1 


sin  0 


where  subscripting  with  r,  6,  or  <f>  denotes  partial  differentiation  with  respect  to  that 
variable.  A  function  u  is  said  to  be  harmonic  in  a  region  D  if  in  that  region  it  satisfies 
Laplace's  equation: 

V2w  =  0. 

In  particular,  if  the  boundary  of  D  is  the  unit  sphere  centered  at  the  origin,  then  using 
the  method  of  separation  of  variables,  the  Legendre  polynomials  arise  naturally  as  part 
of  the  solution. 

In  separation  of  variables,  we  assume  that  the  solution  will  be  of  the  form 


(6) 


u(r,  6,  <t>)=f(r)g(0)h(<p). 


Substituting  (6)  into  (5)  and  dividing  through  by  u,  we  obtain 
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(7) 


/(^)'+7ibr(sin<^+ 


1 


h**  =  0. 


The  first  term,  depending  only  on  r,  must  reduce  to  a  constant  which  we  write  as 
v(v+l).  Substituting  v(v+l)  into  (7)  and  multiplying  through  by  sin20,  we  have 

v(v+l)  sin20  +  -SSI.  {smdge)0+  ■!■  fy,-0. 

Now  we  see  that  the  third  term,  depending  only  on  (/>,  must  also  reduce  to  a  constant; 
call  it  —  in2.  Substituting  —m2  for  this  term  and  simplifying,  we  obtain 

(8)  sin2tf  ge  e  +  sin  6  cos  0  ge  +  [  v(v+l)  sin2fl  -  m2  ]  g  =  0. 

By  the  change  of  variable  x  =  cos  6,  (8)  becomes  via  the  Chain  Rule 


dg 
dx'  dx 


v(v+l)- 


m 


l-x' 


g  =  0. 


For  problems  that  are  radially  symmetric,  the  <f>  dependence  can  be  removed  by  setting 
m  =  0,  leaving 


(9) 


M 


££. 

dx1 


_  -V 


dg 

dx 


+  v(v+l)g  =  0. 


When  v  is  a  positive  integer  n,  (9)  is  recognized  as  the  differential  equation  for  the 
Legendre  polynomials.   [Ref  26:  pp.  210-213] 

An  alternate  approach  to  this  problem  (by  the  method  of  images)  uses  Green's 
functions.  As  motivation,  consider  the  Dirichlet  problem  for  the  unit  circle  in  the  plane, 
which  involves  finding  a  harmonic  function  u(r,  0)  in  the  unit  disk  that  takes  on  pre- 
scribed function  values  f(6)  on  the  boundary  r  =  1.  The  solution  is  given  by 


Lit 


Tin 


f(i)Fr(0-Q')dd\ 


where 


Pr(v)  = 


l-rz 
1— 2r  cos  y+rA 
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is  the  so-called  Poisson  kernel  for  this  problem. 

Similarly,  the  solution  to  the  Dirichlei  problem  for  the  unit  sphere   in  R3 

( V2u  =  0 
\u=f(d,4>)  on  r=[ 

can  be  expressed  as 


u(r,  6,  <f>)  = 


4n 


*  rf(d',<t>')Pr(y)  sinO'dO'df 

n      Jn 


0      •'O 


where 


(10) 


Pfy)  = 


\-r' 


[  1—  2r  cos  y+r 


2-13/2 


and 


cos  y  =  cos  0  cos  0'  +  sin  0  sin  0'  cos(</>  —  0'). 


Writing  (10)  in  terms  of  simpler  functions 


PriY)  - 


1 


1— 2r  cos  y+r 


+  2r  -r- 
2  or 


_    U  1— 2r  cos  y+r 


we  note  the  appearance  of  the  generating  function  for  the  Legendre  polynomials,  with 
x  =  cos  y.  Therefore  the  Legendre  polynomials  are  again  part  of  the  solution.  [Ref.  27: 
pp.  87-89] 

Laplace's   equation   can   also   be   solved   in   a   higher  dimensional   setting.      Let 
x  =  (xx, ...  ,  jc)  e  W.       Note       then       that       we       may       write       x  =  re,       where 


x 


r  =  ||  x  ||  =  yj'xl  +  ...  +  x]    and  %  —  —  =  (£h ... ,  £,)  e  R'  is  a  unit  vector.    A  polynomial 
/i„(x)  is  said  to  be  homogeneous   of  degree /7   if  /t„(ix)  =  i"  hn(x). 

The  Laplacian  operator  acting  on  a  function  u(x)  in  p  dimensions  is  defined  as 


V72  V      ^U 


«=1      dxi 
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As  before,  a  harmonic  function  is  one  which  satisfies  Laplace's  equation  V2pu  =  0.  We 
now  seek  homogeneous  harmonic  polynomials  of  degree  n  in  x  e  W.  It  can  be  shown 
that  there  are  exactly 

2/7  +  o  —  2    (n  +  p  —  3 

linearly  independent  such  solutions,  and  they  can  be  characterized  by  Gegenbaucr  (or 
ultraspherical)  polynomials  C^{tk),  k  =  1,  2, ... ,  N.  The  general  solution  is  given  by 

where  the  spherical  harmonic 

N 


SJti-^A^Cy*    )(fc,iu)), 


k=l 

with  %  suitably  chosen  unit  vectors.     Note  that  if  p  —  3,  then  these  reduce  to  the 
Legendre  polynomials  found  earlier.   [Ref.  6:  pp.  168-183] 

G.     GENETICS  MODELING 

Karlin  and  McGregor  gave  an  interesting  application  of  the  dual  Hahn  polynomials 
to  a  model  in  genetics.  In  this  continuous  time  Markoff  chain  model,  the  dual  Hahn 
polynomials  Rn(  a(x)  ;  y,  S,  /V)  arise  in  the  transition  probability  function  for  the  process. 

The  setting  for  the  model  assumes  N  gametes  of  type  a  or  A  and  gives  a  random 
fertilization  scheme.  The  population  of  either  type  of  gamete  is  affected  by  both  the 
fertilization  process  and  a  mutation  process  whereby  a  gamete  resulting  from  a  mating 
can  mutate  into  the  other  type. 

By  considering  the  conditional  probabilities  for  an  increase  in  population  size  of 
both  gamete  types,  a  stochastic  process  is  defined  which  is  a  classical  birth  and  death 
process.  The  transition  probability  function  for  this  last  process  is  then  cast  in  terms 
of  the  dual  Hahn  polynomials.  The  interested  reader  is  referred  to  [Ref.  14]  for  the  de- 
tails. 
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VI.     BASIC  EXTENSIONS 

A.     BASIC  HYPERGEOMETRIC  SERIES 

In  this  chapter,  we  extend  the  structure  and  some  of  the  results  in  Chapter  IV  to  a 
more  general  level.  This  extension  is  accomplished  by  introducing  a  new  parameter 
called  the  base  to  the  hypergeometric  series.   The  base  q  was  used  by  Heine  in  a  series 

(i-oov)     (iv)(ir')(i-<t)(i-^1)  , 

u         +  d-<,)(i-/)         w(i-V)(i-rt(i^')         " 

where  c  ¥=  0,  —  1,  —2, ...  [Refs.  28,29].  This  series  converges  absolutely  for  |  x  \  <  1  when 
I  q  |  <  1  by  the  Ratio  test.   Since 

1  —  a 
(2)  lim-7— =  a, 

we  see  that  the  series  in  (1)  tends  termwise  to  the  ordinary  hypergeometric  series  as 
<7~>1.  Thus  Heine's  series  is  called  the  basic  hypergeometric  series  or  the  q- 
hypergeometric  series.  [Ref.  12:  p.  3] 

The  q- shifted  factorial  is  the  basic  extension  of  the  shifted  factorial  introduced  in 
Chapter  II  and  is  defined 

fl,  «  =  0 

[(l-a)  {l-aq)  ...  [\-aq      J,    n=\,2,  ...  . 

The  ordinary  shifted  factorial  is  recovered  by  applying  (2)  and  (3)  in  the  limit 

v      ^  :  q^n       i  \ 
lim —  =  (a)„. 

We  also  define 


(4)  («;<?)00  =  I{(l-a/), 

A-=0 
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a  form  wc  will  see  in  later  results.  Since  the  infinite  product  diverges  when  both  a  j=Q 
and  |  q  \  >  1,  we  will  assume  |  q  \  <  1  whenever  (a  ;  q)^  appears  unless  otherwise  stated. 
[Ref.  12:  p.  3] 

Generalizing    the    basic    hypergeometric    series    above,    we    define    the    //>,    basic 
hypergeometric  series  (or  ,<fis  series) 


(5)      A 


au  ...  ,  ar 
\_b\,  •••  i  b. 


",q,x 


_V  {ax\q)n...{ar\q)n  /  Q  \1+w    „ 

Lj    (q;q)n(b];q)n...(bs-q)n     V  {       )    *         )  * 


-l       n=0 


where   (     \—n{n—  l)/2  and  ^  ^  0  when  r>5+l.   We  require  that  the  parameters 
&,, ...  ,  b,  be  such  that  the  denominator  factors  in  each  term  of  the  series  are  nonzero. 


Since 


(<7  m  I  4)n  =  0  .  n=*m+  l,m  +  2t... , 


a  ,0,  series  terminates  if  one  or  more  of  the  numerator  parameters  is  of  the  form  <7_m  for 
m  =  0,  1,2, ...  and  <?  ?t  0.  When  r  =  s  +  1,  the  expression  in  (5)  simplifies  to 


5+1^5 


</>, 


a,,  ...  ,  aJ+1 
ftj,  ...  ,  6j 


;?,* 


-i 

-J        «=0 


(a,  ;^)w...(gf+1  ;<7)„ 
->    ((7 ;  qr)„  (/?,  ;  ?)„  ...  (65 ;  q)n 


Note  that  in  a  basic  hypergeometric  series  X  cn-*">  the  ratio  cn+1/cn  is  a  rational  function 
of?".  [Ref.  12:  p.  4] 

Using  the  q  —  shifted  factorial  (3)  and  (4),  we  can  define  basic  extensions  for  many 
of  the  functions  and  formulas  introduced  in  earlier  chapters.  We  note  that  often  there 
is  more  than  one  way  to  extend  a  result;  examples  will  be  given  shortly. 

The  q-gamma  function   is  defined  by 


(6) 


rq(x)  =  -^f-(i-qy-x,  0<q<l. 


Gosper  showed  that 
(7) 


lim  rxv)  =  r» 


for  0  <  q  <  1  [Ref.  30:  p.  109].    The  structure  of  the  gamma  function  extends  as  well. 
For  instance,  the  formula 
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rxv+i)  = 


i    * 


rxt) 


can  be  reduced  using  (6)  to 


r(.r+i)  =  jcr(*). 


With  (6)  we  can  define  the  q-beta  function 


which  by  (7)  tends  to  B(x,y)  as  q  -*  \~. 

The  q-binomial  coefficient  is  defined  for  integers  n  and  k   by 

(<7 ;  q)n 


Lkl      (q ; 


q)k{q\q)n-.k  ' 


where  k  =  0,  !,...,«.  For  nonintegral  a  and  /?,  we  define 


(/+1;^(^+1;d 


Ak        («? ;  e)oo  (<7a+1 ;  A 
r>+i) 


n/i+i)r((a-j8+i)  • 


The  q-binomial  theorem   is  then 


(a6 ;  ^)n  =  2-,  L.    ^fe  (« ;  ?)*  (£ ;  <?)*-* 


/<=0 


where  «  =  0,  1,2,...  .  [Ref.  12:  p.  20| 

The  next  two  expressions  are  basic  extensions  of  the  Chu-Vandcrmonde  formula 
(Chapter  III,  Section  F.  1)  and  are  both  known  as  the  q-Chu-V cinder monde  formula: 


201 


q     ,  b  cq 

c      '«'  — 


(elb  ;  g)n 
(c  ;  q)n 
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2</M 


—n     i 
q     ,  b 

c 


— : —  b 


[Ref.  12:  p.  11].    These  forms  can  be  shown  to  be  equivalent  by  reversing  the  order  of 
summation. 

Jackson  introduced  the  general  form  of  the  q-integral 


\  f{t)dqt=\  f{t)dqt-\af{t)dq 

Ja  J0  J0 


where 


00 
I    f{t)dqt=a{\-q)Yf{aqn)qn 


The  q  —  integral  defines  the  measure  dqt  which  is  a  natural  object  for 
q  —  defined  functions  [Ref.  12:  p.  19].  All  of  the  functions  and  formulas  developed  above 
play  an  important  role  in  generalizing  the  ordinary  orthogonal  polynomial  classes. 

B.     BASIC  EXTENSIONS  OF  ORTHOGONAL  POLYNOMIALS 

In  this  section,  we  present  the  q  —  analogue(s)  of  selected  classes  from  earlier  chap- 
ters. By  using  formulas  such  as  those  presented  in  the  previous  section  together  with 
methods  based  on  those  outlined  in  the  preceding  chapters,  it  is  possible  to  derive  the 
recurrence  relations,  difference  equations,  and  Rodrigues'  formulas  as  well  as  many 
other  identities  satisfied  by  these  q  —versions.  As  mentioned  in  the  previous  section,  the 
q  —  extension  of  a  function  is  not  necessarily  unique;  however  the  last  two  classes  listed 
are  especially  important. 

1.     Continuous  q  —  Hermite  Polynomials 
a.     Definition 


k=0 

where  jc  =  cos  0. 


Hn(X\q)  =  t    ■        {q)q)n e^-lk)0 
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b.     Orthogonality  Relation 


\*  Hm{  cos  0  |  q)  Hn{  cos  B\q)\  (eW  ;q)J2  JO  =  ^p*    (5 


where  .v  =  cos  6.  [Ref.  12:  p.  188] 

2.     Discrete  <?  — Hermite  Polynomials 
a.     Definition 


[«/2] 


//„(*;</)=£ 


(<7  ;  <7)n 


=o  (?  ;  <i  )k  (q ;  ^)«-2* 


/       n<:     fc(Jfc-l)     /7-2/c 
(-1)    </  ;* 


Z>.     Orthogonality  Relation 


r 


7/wU  ;  <?)  Hn(x  ;  ?)  #(*)  =  qw  (q  \  q)n  ^>m,n 


-aQ 


where  \J/(x)  is  a  step  function  with  jumps 


I  x  |    (■*  <?  ;  <7  )oo  (? ;  <7  )c 


(?  ''i") 


at  the  points  x  =  ±  qJ ,  j  =  0,  1,2,.. 
3.     <7  — Laguerre  Polynomials 
a.     Definition 
For  a  >  —1, 


[Rcf.  12:  p.  193] 


(3) 


CV;<?)  = 


/    a+l        \ 

(g     ;  gj* 
(q ;  <?)« 


i^i 


'7  fl+a+1 

<7 


Z>.     Continuous  Orthogonality 


(y) 


L%{x;q)L™{x;q) 


W( .  • 


.rVv  r(a+l)  f(-a)  (q*     ;$)„ 


(-(l-^v;</)c 


rj-a)  (<7 ;  q)n  ci 


n  ^m,n 


92 


c.     Discrete  Orthogonality 


oo 


1  (a)/       k        \    r  (i)(       k        \ 

L)n  [cq    ;  q)  L\   (cq    J  q) 


/<(3+l) 


/<  =— OO 


=  A 


n      "m,n 


where 


A  = 


(q ;  ?)oo  (-c(i-<7)<7g+1 ;  q)^  (-(i-?)/(cQ ;  ?)c 
(<7a+1 ;  ^)oo  (-cU-<7) ;  ^)oo  (-  —  (i-q) ;  ^)oo 


[Ref.  12:  pp.  194-195] 

4.     Little  #  — Jacobi  Polynomials 
«.     Definition 


p„(x;a,b;q)  =  2(j)] 


q     ,  abq 
aq 


;q,  xq 


b.     Orthogonality  Relation 


2^  pm(qx  \a,b\q)  pn(qx  ;a,b\  q) 


x=0 


(bq  ;  q)x 

(q ;  q)x 


{aq)' 


(q  ;  q)n  (l-abq)  (bq  ;  q)n  (abq1  ;  g)^ 
(a&<7 ;  <7)„  (l-a^2/!+1)  {aq  ;  ?)„  (a? ;  q)c 


(«?)"  ?>m,n 


where  0  <  q  ,  aq  <  1  [Ref.  12:  p.  166] 
5.     Big  q  —Jacobi  Polynomials 
a.     Definition 


Pn{x  ;a,b,c;q)  =  3</>2 


q     ,  abq       ,  x 


aq,  cq 


;<7>  q 
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b.     Orthogonality  Relation 


r 

J cq 


[xja  ;  q)^  {xjc  ;  q)^ 

Pm(x  ;a,b,c;q)  Pn(x  ;  a,  b,  c  ;  q)  — ■ — — — ■ —  d  x 

(*  i  ^)oo  (bxlc  ;  q)^      q 

, ,       (<T,  <l)n  ( 1  ~ah(l)  {bq  ;  q)n  (abq/c  ;  </)„  (") 

=  M : 7—. (— ac)     a\2/ 


(abq  ;  q)n  (\-abq  n+  )  (aq  ;  q)n  {cq  ;  q)n 


(-ac)      q\iJ  dmn 


where 


iV/=         — : — - — : ; —  a„x 


'aq    (xja  I  q)^  {xjc  ;  g)c 

(x ;  ^7)00  (&*/c ;  q)^   "T 


aq(l-q)  (q  ;  ^  (c/a  ;  g)^  fcg/c  ;  q)^  (abq2  ;  g)r 


[Ref.  12:  pp.  167-168] 

6.     #  — Krawtchouk  Polynomials 
a.     Definition 


Kn{x  ;  a,  N;q)  =  3<b2 


—n                 — 
q     ,  x,  -a 

1    n 
1 

-N    A 
q     ,  0 

;q, 

q 

b.     Orthogonality  Relation 


N 


(q    ;  ?)* 


Km(q     \a,N',q)  Kn(q  x  ;  a,  N  ;  q)  (-a) 


x=0 


(<7 ;  q)> 


,     -x    ,    N  -(N+l)   Ur,q)n(i  +  a)(-aqN+l;q)n 

=  (-qa      -q)Na    q   \  2   )  — 

{-a      ;q)n(\+a     q    )  (q       \q)n 


x  (-aq     )    q       tm<n 


[Ref.  12:  p.  185] 

7.     q  — Hahn  Polynomials 
a.     Definition 


Qn(x  ;  a,  b,  N  ;  q)  =  34>2 


n+]        -x 


q     ,  abq       ,  q 

-n      >q>  q 

aq,  q 
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b.     Orthogonality  Relation 


N 


(aq;q)x(bq;q)w_x 

(q ;  i)x  w ;  q)s-. 

(abq2  ;  q)N{aq)~'  (q  ;  ?)„  (l-<%)  (^  ;  ?)„  (abqN+2  ;  ?)„ 


2^  CnC*  ' a' 6- iV  -  ?)  CnC* ;  a>  *.  ^v ;  <?) 


(7  ;  q)s  (abq  ;  ?)„  (l-^2w+1)  (aq  ;  <7)„  (?   ,V  ;  q)n 


-aa)"aO" 


x  (-«?)  q 


■Nn 


Sm„,    m,  w  =  0,  1,  ...  ,  A/" 


[Ref.  12:  p.  165] 

8.     Dual  q  — Hahn  Polynomials 

a.  Definition 

Rn(  n(x)  ;b,c,N;q)  =  3(J>2 

where  fx(x)  =  q~"  4-  cq"~N 

b.  Orthogonality  Relation 


—n        —x         x—  iV 
q     ,  q     ,  cq 

-n  .        >q>  q 

q      ,  bcq 


N 


x=0 


]T  Rm(n(x) ;  b,  c,  N ;  q)  Rn{ji(x) ;  b,  c,  N ;  q) 

[cq~N ;  q)x  ( 1  -cq2x~N)  {bcq  ;  ^)x  (?~'V  ;  ?)x      Nx_  (x\ 

(q;q)x{l-cq      )  (b     q 
(l/c;<7).v         (q;q)n(bq;q)n 


(-bcq)' 


{bq\q)N      (bcq;q)n{q-'W;q)n 


{cq    '  Y  <>m,n  t     m>    «  =  0,  1,  ...  ,  N 


[Ref.  12:  p.  166] 

9.     q  — Racah  Polynomials 
a.     Definition 


(8) 


pn(n(x) ;  a,  b,  c,  d ',  q)  =  4</>3 


"   —  n        1    n+\        —x        t  .T+l 
q     ,   abq       ,   q      ,   cdq 


aq,  bdq,  cq 


;?.  q 


where  /*(*)  =  qr*  +  cdq'A 
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b.     Orthogonality  Relation 

When  aq  ,  bdq  ,  or  cq  =  q~N  for  a  positive  integer  N,   then 


:V 


(9)  2-i PfM*) ;  a'  *'  c'  ^ '  ^ ^(^W  ;a,b,c,d;  q)  w(x  ;  a,  b,c,d;q)  =  hn  Smn 


x=0 


where 


w(x  ;a,b,c,d;  q) 


(10) 


2jc+1 


{cdq  ;  </)x  (ag  ;  q)x  {bdq  ;  <7)x  (07  ;  q)x 1  -  cdq 

{q\q)x  (a~lcdq  ;  q)x  [b~Xcq  ;  q)x  {dq  ;  q)x  {abqf  l  ~  cda 


and 


(11) 


h„  = 


(q  ;  q)n  {\-abq)  (bq  ;  ^  (a^~V  ;  ^  (a6c_1?  ;  </)„  (cafy)" 
(a£<7  ;  </)„  (i-abq2"+])  {aq  ;  <?)„  (^  ;  q)n  {cq  ;  </)„ 


[Ref.  20:  p.  1014].    When  a<?  ,  bdq  ,  or  c<?  is  equal  to  q~N,  the  infinite  products  in  (11) 
reduce  to  finite  products.   Hence  the  orthogonality  relation  (9)  is  valid  for  all  q  provided 
no  zeros  are  introduced  into  denominator  terms  [Ref.  18:   p.  4]. 
10.     Askey- Wilson  Polynomials 
a.     Definition 


pn(x  ;a,b,c,d\  q)  =  a  n  {ab  ;  q)n  {ac  ;  q)n  {ad ;  q)n 


(12) 


x40 


47^3 


n-\  19  -i9 


q     ,  abedq      ,  ae   ,  ae 
ab,  ac,  ad 


;?.  q 


where  x  =  cos  6  [Ref.  18:  p.  3] 

b.     Orthogonality  Relation 
For  —  1  <  a,  b,  c,  d ,  q  <  1 , 


—  J     pm(x  ;  a,  b,c,d\  q)  pn(x  ;  a,  b,  c,  d  |  q) 


w{x  ;  a,  b,  c,  d  \  q)  dx 


\ 


/l- 


"n  "rn,n 
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where 


w 


(x  ;a,b,c,d\  q)  = 


1(1-2(2.^-1)/  +  ^) 


k=0 


h(x,a)  h(x,b)  h(x,c)  h(x,d) 


with 


h(x,a)  =         ( 1  —  2axq    +  a  q    ) 


£=0 


and 


h„  = 


(abcdq2n  ;  q)^  (abcdqn  X  ;  q)n  (qn+]  ;  q)J  [abqn  ;  q) 


-i 


(ac?"  ;  tfL  {adqn  ;  9)^  (fcc?"  ;  q)^  (bdqn  ;  ^  (afy"  ;  q\ 


[Ref.  18:  pp.  11-14] 

We  can  establish  a  formal  connection  between  these  last  two  classes.   If  the 
q  —  Racah  polynomials  are  written 


03 


4^3 


— n       ,  it    n+\        —  x       ,  ,,    x+\ 
q     ,  a  b  q      ,  q     ,  c  d  q 

a'q,  b'd'q,  c'q 


q,  q 


and  the  following  parameter  changes  are  made: 


ac  ,,       bd 

a  =  —     ,    b  =  — 


c  = 


ab_ 

q 


„         u 


(or  likewise  using  any  permutation  of  {  b,  c,  d }  assuming  at  least  one  of  these  is  non- 
zero), as  well  as  the  change  of  variable  q~*  =  ae'e.  then  we  obtain 


403 


i0  -iO 


—n        r     j   n—l  w 

q     ,  abcdq      ,  ae    ,  ae 
ab,  ac,  ad 


q,  q 


with 


=  ae'  +ae  '"  =  2a  cos  0 


/j.(x)  =  q   x  +  qx     c'd 
iO  ,      -iff 
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That  is,  the  q  —  Racah  and  Askey-Wilson  polynomials  arc  virtually  the  same,  differing 
only  in  their  parameters,  normalization,  and  variable. 

Because  the  q—  Racah  and  Askey-Wilson  polynomials  are  essentially  the 
same,  the  names  are  often  used  interchangeably  in  the  literature.  We  have  presented 
these  polynomials  as  distinct  classes  in  order  to  emphasize  the  continuous  and  discrete 
natures  as  expressed  in  the  orthogonality  relations. 

As  stated  at  the  beginning  of  Chapter  IV,  the  "classical  orthogonal 
polynomials"  are  defined  to  be  those  which  are  special  or  limiting  cases  of  the  Askey- 
Wilson  (21)  or  (/—Racah  polynomials  (18)  [Rcf.  31:  p.  57].  We  can  now  make  this 
statement  a  bit  more  precise.    Letting 

a  =  q      ,    b  =  qy    ,    c  =  q      ,    a  =  q 

in  (9)  and  (10)  and  taking  the  limit  as  q->\,  the  ordinary  Racah  polynomials 
R„(  X(x) ;  a,  /?,  y,  <5)  and  their  weight  function  as  defined  in  Formulas  (30)  and  (32), 
Chapter  IV,  Section  C.3  are  retrieved.  A  similar  limiting  process  will  recover  any  ordi- 
nary orthogonal  polynomial  class  from  its  (/—extension.  Moreover,  any  (/—extension 
previously  discussed  is  a  special  or  limiting  case  of  the  q  —Racah.  For  example,  letting 
cq  =  q~N  and  d=0  produces  the  q —Halm,  etc.  In  this  way,  we  see  that  the  self-dual 
q  —Racah  polynomials  encompass  all  the  previous  classes.  They  also  satisfy  three-term 
recurrence  relations,  second  order  difference  equations,  and  Rodrigues'  formulas  with 
respect  to  (/—divided  difference  operators.  The  interested  reader  is  referred  to  [Refs.  18, 
12,  20]  for  details  about  these  very  rich  classes. 

C.     CONCLUDING  REMARKS 

In  Chapter  V,  we  presented  a  few  of  the  traditional  applications  in  which  orthogonal 
polynomials  have  played  an  important  part.  More  recently,  an  enriching  interplay  has 
developed  between  the  theory  of  orthogonal  polynomials  and  other  mathematical  and 
mathematically-related  areas. 

Efficient  computational  methods  have  been  devised  for  determining  the  many  useful 
quantities  associated  with  orthogonal  polynomials  (such  as  their  zeros,  recurrence  coef- 
ficients, etc.)  [Rcf.  32:  pp.  181-216).  Various  classes  of  orthogonal  polynomials  have 
also  played  a  role  in  digital  signal  processing  [Rcf.  32:  pp.  1 15-133],  quantum  mechanics 
[Rcf.  32:  pp.  217-228],  and  birth/death  processes  [Ref.  32:  pp.  229-255].  Advances  in  the 
field  of  combinatorics  and  graph  theory  have  allowed  new  geometric  interpretations  of 
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orthogonal  polynomial  identities,  some  of  which  have  very  important  consequences  for 
"association  schemes"  and  the  designs  of  codes  [Refs.  33,  16,  34,  32  :  pp.  25-53,  35]. 

Physicists  have  introduced  various  versions  of  "diagrammatic  methods":  ways  of 
understanding  orthogonal  polynomials  through  3-j  and  6-j  symbols,  their  generaliza- 
tions, and  accompanying  identities  by  formally  associating  them  with  pictorial  schemat- 
ics representing  forces  of  physical  systems  that  conserve  angular  momentum  [Ref.  36]. 
Powerful  new  techniques  involving  "quantum  groups"  have  been  used  to  generate  new 
identities  for  some  classes  [Ref.  32:  pp.  257-292].  Further  investigation  into  the 
electrostatics  problem  discussed  in  the  text  has  led  to  the  formation  of  the  famed 
"Selberg  beta  integral"  and  its  generalizations.  It  has  yielded  to  analysis  via  the  study 
of  the  "root  systems"  of  Lie  algebras,  and  has  found  applications  ranging  from  statistical 
mechanics  to  computer  algorithm  complexity  [Rcfs.  30:  pp.  48-52,  32  :  pp.  311-318,  37]. 

Finally,  research  into  the  general  structure  of  q  —series  has  led  to  many  surprising 
connections,  and  is  intimately  related  to  the  many  remarkable  and  powerful  number- 
theoretic  formulas  discovered  by  S.  Ramanujan,  the  famed  Indian  mathematical  prodigy 
[Refs.  30:  pp.  87-93,  38:  pp.  55-66,  39]. 

This  utility  and  promise  of  future  applications  provide  ample  motivation  and  justi- 
fication for  continued  study  of  the  intrinsic  structure  of  orthogonal  polynomials. 
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